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We consider in detail the situation of applying a time dependent external magnetic field to a *'Rb atomic Bose- 
Einstein condensate held in a harmonic trap, in order to adiabatically sweep the interatomic interactions across a 
Feshbach resonance to produce diatomic molecules. To this end, we introduce a minimal two-body Hamiltonian 
depending on just five measurable parameters of a Feshbach resonance, which accurately determines all low 
energy binary scattering observables, in particular, the molecular conversion efficiency of just two atoms. Based 
on this description of the microscopic collision phenomena, we use the many-body theory of T. Kohler and 
K. Burnett [Phys. Rev. A 65, 033601 (2002)] to study the efficiency of the association of molecules in a ^'Rb 
Bose-Einstein condensate during a linear passage of the magnetic field strength across the 100 mT Feshbach 
resonance. We explore different, experimentally accessible, parameter regimes, and compare the predictions of 
Landau-Zener, configuration interaction, and two level mean field calculations with those of the microscopic 
many-body approach. Our comparative studies reveal a remarkable insensitivity of the molecular conversion 
efficiency with respect to both the details of the microscopic binary collision physics and the coherent nature of 
the Bose-Einstein condensed gas, provided that the magnetic field strength is varied linearly. We provide the 
reasons for this universality of the molecular production achieved by linear ramps of the magnetic field strength, 
and identify the Landau-Zener coefficient determined by F.H. Mies et al. [Phys. Rev. A 61, 022721 (2000)] as 
the main parameter that controls the efficiency. 

PACS numbers: 03.75.Kk, 34.5().-s, 05.30.-d 



I. INTRODUCTION 



The prospect of achieving quantum degenerate molecular 
gases has attracted considerable attention for some time now 
Ql]. Such an accomplishment may open new avenues for re- 
search, for instance, bright sources of molecules for cold col- 
lision studies 01, precise molecular spectroscopy, elucidat- 
ing the nature of a possible BCS-BEC crossover in Fermi 
gases 0] and, possibly, the exploitation of dipole-dipole in- 
teractions ||4|. It has been clear from the outset, however, 
that laser cooling techniques, essential for the production 
of Bose-Einstein condensates and degenerate Fermi gases of 
atoms, are difficult to apply in the case of molecules due 
to their typically complicated rovibrational energy spectrum. 
The association of ultracold atoms into diatomic molecules, 
which may also be quantum degenerate, therefore seems a 
very promising route. The molecular conversion can be 
achieved by photoassociation |5] or by application of time 
dependent magnetic fields in the vicinity of Feshbach reso- 
nances 1 6]. The Feshbach resonance technique has recently 
been exploited with great success to produce large, ultra- 
cold assemblies of diatomic molecules, usin^ as__a source 
both atomic Bose-Einstein condensates ff, 8, 9, 10,1^ and 
quantum degenerate two component gases of Fermionic atoms 

The experiments reported in Refs. jgl lTollTlll achieve con- 
version of alkali atoms in dilute Bose-Einstein condensates to 
diatomic molecules by an adiabatic sweep of the strength of 
a homogeneous magnetic field across a Feshbach resonance 
from negative to positive scattering lengths. The observations 



indicate molecular fractions much smaller than the ideal limit 
of half the number of initial condensate atoms, even in the 
limit of perfect adiabaticity. This may be due to limitations in 
the initial production of molecules, or their long term stability 
in the presence of the surrounding gas. 

The long term stability of the highest excited vibrational 
molecular bound state, produced by the adiabatic associa- 
tion technique, may be limited due to collisional deexcitation. 
Both theoretically and experimentally, very little is known 
about the associated rate constants for highly excited diatomic 
molecules composed of alkali atoms. To our knowledge, the 
only available exact calculations have been performed for the 
deexcitation of tightly bound Na2 molecules upon collision 
with a Na atom |24]. From an experimental viewpoint, the 
comparatively high relative momenta of the products of colli- 
sional deexcitation have so far prevented their detection. Con- 
clusive experimental studies of the loss rates and the underly- 
ing microscopic processes are therefore difficult to achieve. 

In such an uncertain situation it is important to first under- 
stand the production of the diatomic molecules. In this pa- 
per we consider the situation of adiabatically sweeping the 
magnetic field strength across a Feshbach resonance to form 
molecules from an initially Bose-Einstein condensed dilute 
gas of atoms. The calculations presented in this pap er con- 
sider the 100 mT Feshbach resonance of ^^Rb [251. The 
underlying concepts can be extended to arbitrary species of 
Bosonic atoms, while the two-body physics is also applica- 
ble to any pair of atoms interacting via s waves, including 
Fermionic species. 

The structure of the paper is as follows: SectionHIIdescribes 
the binary physics of the association of ultracold atoms to a 
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diatomic molecule, and Section HiH analyses the many -body 
aspects of this process in a Bose-Einstein condensate. Section 
IIVI summarises and presents our conclusions, and there is an 
appendix, which expands upon the explicit description of the 
binary collision physics specific to this work. 

Section ini presents a two channel description of the mag- 
netic field tunable resonance enhanced binary scattering and 
its relationship to the properties of the highest excited vibra- 
tional bound state. We develop the concept of a two-body 
Hamiltonian that accurately describes the relevant physics in 
terms of a minimal set of five parameters, which can usually 
be deduced from measurable properties of a Feshbach reso- 
nance. We then focus on the near resonance universal proper- 
ties of the highest excited vibrational bound state, and identify 
the smooth transition between free and bound atoms in an ide- 
ally adiabatic passage across a Feshbach resonance. We then 
consider the dynamics of the association of just two atoms 
during a linear ramp of the magnetic field strength across a 
Feshbach resonance, and reveal the dependence of molecular 
conversion on the physical parameters. The last subsection 
is concerned with the determination of dissociation energy 
spectra, as dissociation of molecules by ramps of the mag- 
netic field strength is frequently a necessary precursor to their 
detection. We show that both the molecular conversion efli- 
ciency and the dissociation spectra are determined by the same 
physical parameters of a Feshbach resonance, and are subject 
to a remarkable insensitivity with respect to the details of the 
binary collision dynamics, provided that the magnetic field 
strength is varied linearly. 

Section |ni| introduces the relevant elements of the micro- 
scopic quantum dynamics approach |26], which properly ac- 
counts for both the microscopic binary collision physics and 
the macroscopic coherent nature of an inhomogeneous atomic 
Bose-Einstein condensate. Using this approach, we study the 
many-body aspects of the molecular production during a lin- 
ear adiabatic passage of the magnetic field strength across a 
Feshbach resonance. From the first order microscopic quan- 
tum dynamics approach, by applying the Markov approxi- 
mation, we derive the commonly used two level mean field 
model iHllllllllSlllllllllllllll. We discuss the 
deficits of two level models with respect to the description 
of the intermediate dynamics during passage of the magnetic 
field strength across a Feshbach resonance. Our comparative 
studies show, however, that the many-body approaches pre- 
dict virtually the same final molecular production at all levels 
of approximation considered in this paper, provided that the 
ramp of the magnetic field strength across the Feshbach res- 
onance is linear We provide the reasons for this universal- 
ity of the molecular conversion in linear ramps of the mag- 
netic field strength with respect to the details of the underly- 
ing microscopic binary collision dynamics and to the coher- 
ent nature of the Bose-Einstein condensed gas. We then show 
that the molecular production efficiency is determined by the 
same physical parameters that we have previously identified 
in the associated two-body problem. Our findings strongly in- 
dicate that measurements of molecular production efficiencies 
as well as dissociation spectra obtained from linear ramps of 
the magnetic field strength are largely inconclusive with re- 



spect to the details of the underlying binary collision physics. 



II. ADIABATIC ASSOCIATION OF TWO ATOMS 

We consider a configuration of two atoms exposed to a 
homogeneous magnetic field whose strength can be varied. 
The concept underlying the adiabatic association of diatomic 
molecules in ultracold gases can be understood solely on the 
basis of binary collision physics. The key feature of this ex- 
perimental technique is the adiabatic transfer of the zero en- 
ergy binary scattering state into the highest excited diatomic 
vibrational bound state. In this section we shall study the bind- 
ing energies of the diatomic molecules that determine the po- 
sitions of the Feshbach resonances. We shall then show how 
the resonance enhanced interatomic collisions can be accu- 
rately described in terms of a minimal set of five quantities 
that can be determined from current experiments. Specifi- 
cally, we consider the association of two asymptotically free 
ultracold ^^Rb atoms with a total angular momentum quan- 
tum number of F - 1 and an orientation quantum number 
of nif - -Hi, at magnetic field strengths close to the broadest 
Feshbach resonance at about 100 mT. We shall also study the 
dissociation of the molecules, which plays an important role 
in the direct detection of ultracold molecular gases. 



A. Feshbach resonances and vibrational bound states in *''Rb 

Throughout this paper, we will denote the open scattering 
channel of two asymptotically free atoms in the (F = 1, m/r = 
-Hi) electronic ground state as the open channel with an as- 
sociated reference potential (the background scattering poten- 
tial) Vbg{r). The dissociation threshold of Vbgir) is determined 
by the internal energy of the noninteracting atoms, i.e. twice 
the energy corresponding to the (F - l,mf = -Hi) hyper- 
fine state. When the atoms are exposed to an external ho- 
mogeneous magnetic field the nif degeneracy of the atomic 
hyperfine levels is removed by the Zeeman effect. As a con- 
sequence, the potentials associated with the different asymp- 
totic binary scattering channels are shifted with respect to each 
other Although in general the interchannel coupling is weak, 
in the vicinity of certain magnetic field strengths the open 
channel can be strongly coupled to closed channels. This 
strong coupling leads to singularities of the s wave scattering 
length known as Feshbach resonances. Figuren(a) shows the 
theoretically predicted s wave scattering length of two collid- 
ing ^^Rb atoms in the electronic ground state that are exposed 
to a homogeneous magnetic field of strength B. The theoret- 
ical calculations use five coupled equations for one open and 
four closed channels to describe the i-wave collision of two 
(F = l,mf = -Hi) atoms f3^. We use standard methods to 
solve these equations for either scattering or bound states. 

The Feshbach resonances are related to the binding energy 
of the highest excited vibrational state in a simple way; the 
singularities of the scattering length in Fig.^Ca) exactly match 
those magnetic field strengths that correspond to the zeros 
of the binding energy with respect to the threshold energy 
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for dissociation into two asymptotically free atoms. Figure 
[0 (b) gives an overview of the magnetic field dependence of 
the binding energies of s-wave symmetry molecular states of 



two Rb atoms below the dissociation threshold energy of the 
open channel. 
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FIG. 1: Magnetic field strength dependence of the s wave scattering 
length (a) and the energies of the highest excited i-wave vibrational 
molecular bound states (b) of two '^'Rb atoms. At each magnetic field 
strength B the zero of energy is set at the threshold for dissociation 
into two asymptotically free atoms in the {F = l,mf = +1) hyper- 
fine state. The lines in (b) at about -0.02 and -0.63 GHz parallel to 
the £ = axis represent the energies of the last j-wave bound states 
of the open background channel. These states have the same mag- 
netic moment as the separated atoms. The slanted lines in (b) rep- 
resent closed channel molecular states that have different magnetic 
moments from the separated atoms. The weak avoided crossings at 
the intersection of parallel and slanted lines are due to interactions 
between the background and closed channels, as noted in Ref. floh . 
When the closed channel molecular states cross threshold at £ = 0, 
the zeros of the binding energy in (b) correspond to the positions of 
the Feshbach resonances in (a). Our calculated resonance positions 
are within one per cent of the measured positions L25il . 



B. Two channel energy states 

The scattering lengths and binding energies in Fig.[2have 
been obtained from exact solutions of the multichannel two- 
body Schrodinger equation as described in Subsection 111 Al 
These calculations were performed with a realistic potential 



matrix that accurately describes the bound and free molecu- 
lar states over a wide range of energies and magnetic field 
strengths. Based on these exact considerations, Fig.^reveals 
that the binding energy of the highest excited vibrational state 
determines the singularities of the scattering length and, in 
turn, all resonance enhanced low energy scattering proper- 
ties of two atoms. Collisions in ultracold gases involve only 
a quite limited range of energies, and the adiabatic associa- 
tion of molecules takes place at magnetic field strengths in 
the close vicinity of a particular Feshbach resonance. We 
shall therefore restrict our analysis to an appropriately smaller 
range of energies and magnetic field strengths around the 100 
mT Feshbach resonance. 



1. Background scattering 

We shall first show how the low energy background scat- 
tering can be accurately described in terms of experimen- 
tally known physical quantities. At magnetic field strengths 
asymptotically far from the resonance the interchannel cou- 
pling is weak and the highest excited multichannel vibra- 
tional bound state can be determined directly from a single 
channel description with the background scattering potential 
Vbg(r). Considerations beyond the scope of this paper iIjtIi 
show that the coiTesponding binding energy is determined, to 
an excellent approximation, by the long range asymptotic be- 
haviour of Vbgir) and the background scattering length abg, 
i.e. the scattering length associated with the background scat- 
tering potential Vbg{r). Neglecting retardation phenomena, 
at large interatomic separations Vbgir) has the universal form 
^bg('") ~ -Ce/r^, where Ce is the van der Waals dispersion 
coeflicient. The low energy background scattering is deter- 
mined by the same parameters Obg and Ce iSTil . This univer- 
sality is due to the fact that at typical ultracold collision ener- 
gies the de Broglie wavelengths are much larger than the van 
der Waals length 
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(1) 



Here m is the atomic mass. As the van der Waals length is 
the characteristic length scale set by the long range tail of 
the background scattering potential, the details of the poten- 
tial Vbgir) are not resolved in the collisions. At zero collision 
energy the background scattering length Obg incorporates all 
the unresolved details of Vbgir) into a single length scale. The 
van der Waals length determines the first correction at finite 
collision energies, which accounts for the long range asymp- 
totic behaviour of Vbgir). We note that any model of the back- 
ground scattering potential will recover the binding energy of 
the highest excited vibrational state and all low energy scat- 
tering properties of the exact potential Vbgir) to an excellent 
approximation, if it properly accounts for the parameters Ce 
and Obg. We provide an appropriate minimal background scat- 
tering potential in the appendix. 
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2. Two channel Schrddinger equation 

We assume in the following that all the bound and free en- 
ergy states associated with the background scattering potential 
have been determined, and on this basis derive the resonance 
enhanced collision properties of two atoms. In the vicinity of 
a Feshbach resonance, the strong coupling between the open 
channel and other asymptotic scattering channels originates 
from the near degeneracy of the magnetic field dependent en- 
ergy EresiB) of a closed channel vibrational state (a Feshbach 
resonance level) 0re.s('") with the dissociation threshold energy 
of the open channel. Consequently, the resonance enhanced 
collision physics of two **^Rb atoms can be accurately de- 
scribed by the general form of a two-channel Hamiltonian ma- 
trix of the relative motion of the atoms: 



H2B - 



+ ybg(r) 

W(r) 



W{r) 
-fv2 + yd(B,r) 



(2) 



Here m is the atomic mass, r is the distance between the 
atoms, W{r) determines the strength of the coupling between 
the channels, and the closed channel potential Vci{B, r) sup- 
ports the resonance state: 



v2 + y,,(B,r) 



(3) 



In the following the resonance state (pi-es{r) is normalised to 
unity. In Eq. Q, as elsewhere in this paper, we have cho- 
sen the zero of energy as the dissociation threshold of the 
open channel, i.e. V'bgC'") 0. The dissociation threshold of 

r— >oo 

VcliB, r) is determined accordingly by the energy of two non- 
interacting atoms in the closed channel that is strongly cou- 
pled to the open channel. The relative Zeeman energy shift 
between the channels can be tuned by varying the magnetic 
field strength. 

The bound and free energy states of the general Hamilto- 
nian matrix in Eq. (|5Jl relate the remaining potentials W{r) and 
Vci{B, r) to a minimal set of measurable properties of a Fesh- 
bach resonance. The two channel states that we shall consider 
in the following are of the general form |bg)0''°(r) + |cl)0'^'(r), 
where |bg) and |cl) denote the internal states of an atom pair 
in the open channel and the closed channel strongly coupled 
to it, respectively. The two components of the energy states 
are solutions of the stationary coupled Schrodinger equations 



V2 + Vt,(r) 

m 



(r) + W(r)0"'(r) = ^^''^^(r), (4) 



— V2 + y,i(B,r) 
m 



Ar) = E<p'\r). (5) 



3. Continuum states 

Bound and continuum energy states are distinguished by 
their energies and by their asymptotic behaviour at large in- 
teratomic distances. In the scattering continuum above the 
dissociation threshold all energies are in the spectrum of the 



Hamiltonian in Eq. (|2} and can be associated with the momen- 
tum p of the relative motion of two asymptotically noninter- 
acting atoms in the open channel through E = p^/m. Due to 
the continuous scattering angles between the atoms at a defi- 
nite collision energy, the scattering energy states are infinitely 
degenerate. In the following, we will choose their wave func- 
tions to behave at large interatomic distances like: 



1 



(6) 



This long range asymptotic behaviour corresponds to an in- 
coming plane wave and an outgoing spherical wave in the 
open channel. Here and throughout this paper we will as- 
sume the plane wave momentum states to be normalised as 
exp(/p ■ rlK)l{2nKfl^. The function /()?, p) in Eq. (|6|l is the 
scattering amplitude, which depends on p - 'sfmE and on the 
scattering angle & between the momentum p of the relative 
motion of the asymptotically noninteracting incoming atoms 
and their final relative position r. The closed channel com- 
ponent 0p (r) of the wave function vanishes at asymptotically 
large distances between the colliding atoms. We shall also in- 
troduce the energy states <p'^\r) of the background scattering 
that satisfy the Schrodinger equation 



ft2 



-v^ + ybg(r) 



(r) 



(r). 



with the long range asymptotic behaviour: 
1 



e'^'^"+f,s(i^,p) 



Jpr/n 



(7) 



(8) 



The coupled Schrodinger equations (0} and Q can be ex 
pressed in terms of the energy dependent Green's functions: 



Ghgiz) 



GAB,z)^ 



— v^ + ybj 

m 



z-\--V^ + VAB) 
m 



(9) 



(10) 



Here z is a complex parameter with the dimension of an en- 
ergy. The coupled Schrodinger equations then read: 



+ Gbg(£ + iO)Wf^, 



Ga{B,E)Wf/ 



(11) 
(12) 



The argument "z - E + /O" of the Green's function Gbg(z) in- 
dicates that the physical energy E - p^/mis approached from 
the upper half of the complex plane. This choice of the energy 
argument ensures that the scattering wave function 0p^(r) has 
the long range asymptotic form of Eq. (|6j, in accordance with 
the asymptotic behaviour of the Green's function at large in- 
teratomic distances: 



Gbg(E + iO, r, r') 



-(Inn) 



3/2 _ 



[cll-\r')] . (13) 



Here 0p '(r') = |^(/>*!y(r')] is the incoming continuum energy 
state associated with the background scattering [38,1 . and p = 
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i^yfmE^rlr can be interpreted as the asymptotic momentum 
associated with the relative motion of the scattered atoms. 

As the resonance state (p^sir) fulfils the Schrodinger equa- 
tion (|3|l, according to Eq. ( I10> the Green's function Gc\{B,z) 
has a singularity at z = E^^^iB), i.e. 



<</'res|Gcl(B,z)l<^res) = 



1 



(14) 



At magnetic field strengths in the vicinity of a Feshbach 
resonance E^^AB) is nearly degenerate with the dissociation 
threshold energy of the open channel. Furthermore, the ki- 
netic energies E - p^/m in ultracold collisions are small in 
comparison with the typical spacing between molecular vi- 
brational bound states. As a consequence, the denominator 
in Eq. (I14> becomes sufficiently small for the Green's func- 
tion Gc\{B,E) in Eq. (I12> to be excellently approximated by 
its resonance state component L39.I : 



Gd(B,£)^ |<^re,s> 



1 



E - Erss(B) 



(15) 



Inserting this pole approximation of the Green's function into 
Eq. (I12t determines the functional form of the closed channel 
component of the scattering wave function to be 



The wave function (^p^(r) is then determined by eliminating 
ipp on the right hand side of Eq. il 1> in terms of Eq. (I16> 
which gives 



|0res)A(B,£). 



(16) 



l-^p") = l<*> + Gbg(£ + iO)W\cf>r,,}A(B,E). 
The as yet unknown amplitude 



A{B,E) : 



(17) 



(18) 



E - £'res(B) 

can be determined straightforwardly by multiplying Eq. ( I17> 
by {(piealW from the left. This yields, after a short calculation: 



The scattering length a is determined by the long range 
asymptotic form of the scattering wave function 0p°(r) in 
Eq. (|6j through 

f(§,p) ~-a. (21) 



Because the incoming plane wave is isotropic at zero energy, 
this limit is obviously independent of the scattering angle. In 
an analogous way, the background scattering length a^o is re- 
lated to the asymptotic form of the scattering wave function 
0p^'(r) in Eq. (|8j, in the limit of zero energy. Inserting the 
known asymptotic form of the Green's function G\,giE + /O) 
at large interatomic distances in Eq. M3\ and the amplitude 
in Eq. into Eq. a short calculation determines the 
scattering length to be: 



£,es(B) + <0i-es|WGbg(O)H^|0,,s) 



(22) 



The energy E^^^iB) of the closed channel state is determined 
by the Zeeman energy shift between the asymptotic open 
channel and the closed channel strongly coupled to it. Within 
the limited range of magnetic field strengths that we shall con- 
sider in this paper, the Zeeman eff'ect is approximately linear 
in the magnetic field strength B. We shall denote by Z?res the 
magnetic field strength at which E^^^{B) crosses the dissoci- 
ation threshold of the open channel, i.e. £'ies(Bres) = 0. An 
expansion of E^esiB) about B = Bres yields: 



ErJB) 



dEr, 



dB 



-(Bks) 



(B - Bres). 



(23) 



Equations J22l i and j23> then determine the magnetic field de- 
pendence of the scattering length by the well known formula 



aiB) = flbg 



(AB) 



B-Bo 



(24) 



A{B,E) = 



E - Eres{B) ■ 



s|WGbg(£' + iQ)W\^ res/ 



(19) 



Once all the energy states associated with the single chan- 
nel potential ybg('") and the resonance state 0res('") are known, 
Eqs. ( I16> , (O and jl9> establish the complete solution of the 
coupled Schrodinger equations (II 1> and il2\ in the pole ap- 
proximation. Under the assumption that the configuration of 
two atoms in the closed channels is restricted to the resonance 
state 01-es, the pole approximation becomes exact. This as- 
sumption implies the replacement 



V2 + VdB) 

m 



l0i-es)£i-es(B)<0 n 



(20) 



in the two channel Hamiltonian in Eq. 0. Equations ( I16> , 
and il9\ then determine the exact scattering energy states 
of the resulting restricted two channel Hamiltonian. 



where 



(AB) 



=,siwi0[;') 



(25) 



is termed the resonance width and 



B() - B„ 



s res/ 



(26) 



is the measurable resonance position, i.e. the magnetic field 
strength at which the scattering length has a singularity. We 
note that according to Eq. i26i the resonance position Bq is 
shifted with respect to the magnetic field strength Bres at which 
the energy of the resonance state becomes degenerate with the 
dissociation threshold energy of the open channel. 
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4. Bound states 

The molecular bound states vanish at asymptotically large 
interatomic distances, and their energies are below the disso- 
ciation threshold of the open channel. Similarly to Eqs. M W 
and M2\ . the coupled Schrodinger equations and (jSj can 
be expressed in terms of the coupled integral equations 



(27) 
(28) 



which incorporate the long range asymptotic behaviour of the 
components, 0]^° and of the bound state. Here E\, is the 
binding energy. A short calculation verifies that the pole ap- 
proximation in Eq. M5\ leads to the normalised solutions 



1 



Gbg(£b)W0„ 



(29) 



with the normalisation constant 



/^b = ^ 1 + <<^res|W [Gbg(£b)]' W|<^,es), (30) 



whenever the binding energy fulfils the condition: 

£b = EM + <<^res|H'Gbg(£b)W|<^res). 



(31) 



We note that Eq. (13 U recovers Eq. \26\ when the binding en- 
ergy and the magnetic field strength are inserted as Zsb = and 
B - Bo, respectively, i.e. the binding energy, indeed, vanishes 
at the position of the resonance. 



5. Minimal two channel Hamiltonian 

Both the resonance width in Eq. (125 > and the shift in 
Eq. (I26> depend only on the product W(r)<piss(r). Conse- 
quently, for a minimal description of the resonance enhanced 
scattering we only need to specify W(r)0i-es(r) in terms of two 
parameters to recover the magnetic field dependence of the 
scattering length in Eq. (I24> and its relationship with the bind- 
ing energy of the highest excited vibrational state. A deriva- 
tion beyond the scope of this paper shows that, once the width 
of the resonance is known experimentally, a determination of 
the resonance shift does not require a full solution of the cou- 
pled channel two-body Schrodinger equation with a realistic 
potential matrix. It turns out that Eq. ( I26t is excellently ap- 
proximated by i4(^l 



Bo - Bres = (AB) 



a \ a I 



(32) 



Here a is the mean scattering length of the background scat- 
tering potential V'bg('"), which is related to the van der Waals 
length in Eq. (^ and Euler's F function by i4lll 



- 1 r(3/4) , 



(33) 



This approximation is consistent with the treatment of the 
background scattering in terms of just the parameters Obg and 
Zvdw in lllB II 

In addition to the parameters C5 and cibg of the background 
scattering potential ybg('"), a minimal description of the reso- 
nance enhanced scattering in a two channel Hamiltonian thus 
requires us to account for the width (AB) in Eq. \25\ and the 
shift Z?o - Bres in Eqs. M6\ and \'i2l . in terms of the inter- 
channel coupling VK(r)0res('"), while the slope [^^^(^les)] of 
the resonance in Eq. ( I23> determines the component of the 
Hamiltonian in the closed channel that is strongly coupled to 
the open channel. For the 100 mT Feshbach resonance of 
**^Rb, four of the five parameters of the two channel Hamil- 
tonian, i.e. Ce = 4660 a.u. |42] (1 a.u. = 0.095734 yj nm*'), 
flbg = 100 flBohr i4i ((JBohr = 0.052918 nm), (AB) = 0.02 mT 
and B() - Bres = -0.006371 mT, can be either directly deduced 
from experiments i44ll or from Eq. \32\ . The only parameter 
that is not easily accessible is the slope of the resonance. We 
have obtained \ [^(Bres)] = 38 MHz/mT from the binding 
energies in Fig. ^ In the appendix we provide the explicit 
form of the minimal two channel Hamiltonian in the pole ap- 
proximation that we apply in the following to determine the 
dynamics of the adiabatic association of molecules. Figure|2 
shows the magnetic field dependence of the binding energies 
as obtained from this minimal Hamiltonian. The low energy 
scattering properties of two asymptotically free ^^Rb atoms 
in the open channel and the properties of the highest excited 
vibrational bound state are insensitive with respect to the de- 
tails of the implementation of the five parameter two channel 
Hamiltonian. 




59 99.5 100 100.5 101 101.5 102 
B[mT] 



FIG. 2: The magnetic field dependence of the binding energy of the 
vibrational bound states of the two channel Hamiltonian in the ap- 
pendix (solid curves). On the low field side of the resonance a new 
bound state 0b emerges at the resonance position Bo = 100.74 mT 
1 44] whose binding energy we have denoted by Eb(B). At magnetic 
field strengths asymptotically far from the resonance, the highest ex- 
cited two channel vibrational bound state becomes identical with the 
highest excited single channel vibrational bound state <p-i(r) that is 
associated with Vbg(r) and whose binding energy is denoted by £_! 
(dotted horizontal line). The dashed line indicates the energy £'|.es(-6) 
of the resonance state. 
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C. Universal properties of near resonant bound states 

We shall now focus on the properties of the highest excited 
vibrational bound state 0|, on the low field side of the Fes- 
hbach resonance whose emergence at the resonance position 
Bo = 100.74 mT causes the singularity of the scattering length 
(see Fig.O- The vibrational state ^t, is determined by its com- 
ponents in the open channel and in the closed channel strongly 
coupled to it, as given by Eq. j29l l. while the binding energy 
is determined by Eq. OH . We note that the closed channel 



component <p'^ {r) in Eq. i29\ has the functional form of the 
resonance state 0res('"), which we have normalised to unity. 
Figure|3shows the population 



47T 



r 

Jo 



r^dr 



1 



(34) 



of the open channel component of the highest excited vi- 
brational bound state 0b as a function of the magnetic field 
strength. 




FIG. 3: The population of the open channel component of the highest 
excited vibrational bound state 0(, as a function of the magnetic field 
strength on the low field side of the Feshbach resonance. The popu- 
lation was determined from the minimal two channel Hamiltonian in 
the appendix. 

At magnetic field strengths asymptotically far from the res- 
onance. Fig. |3l shows that 0b is transferred entirely into the 
bound state (p-i{r) of Vbo{r), i.e. the highest excited vibrational 
bound state in the absence of interchannel coupling. Figure|3 
also reveals that (p\, is dominated by its component in the open 
channel also in a small region of magnetic field strengths in 
the close vicinity of the Feshbach resonance. 



involve the asymptotic behaviour of Eq. (13 H in the limit E\, 
0. Inserting the resolvent identity 138,1 



Gbg(£b) = Gbg(O) - £bGbg(0)Gbg(£b) 
as well as Eqs. j23l l and ( I26I I into Eq. i31\ yields: 



(35) 



-(fires) 



{B - Bo) 



dB 

- £b<0res|WGbg(O)Gbg(£b)W|0res). 



(36) 



In accordance with Eq. the Green's functions Gbg(O) and 
Gbg(£'b) in Eq- (I36> can be decomposed into the complete or- 
thogonal set of bound and continuum energy states associated 
with the background scattering potential ybe('")- These de- 
compositions yield: 



dE,-, 



dB 



■(Bies) 



(B - Bo) 



s|W|<') 



_£:(£,_£) ^-EAE^-E,) 

m \ " in ) ' 



(37) 



Here the sum includes the indices v = -1, -2, -3, ... of all 
vibrational bound states 0,(0 associated with the potential 
Vb^ir). In the limit of vanishing binding energy E^, the mo- 
mentum integral on the right hand side of Eq. ( I37t is singular 
at p = 0. As a consequence, the slowly varying matrix ele- 
ment (^reslW'l^p^*) can be evaluated at /? = 0, while the sum 
over the vibrational bound states can be neglected. Using the 
general relationship between the matrix element (^reslW^I^J,^') 
and the resonance width in Eq. \25\ . the asymptotic form of 
the remaining momentum integral is then given by: 



,s|W|0«) 



m \ " m ) 



An 



dEri 



f 

Jo 



p dp 



siwi0{;') 



m \ m J 



dB 



-(Bies) 



(Ai?)flbgJ^ 



'|£bl 



(38) 



With this evaluation of the integral, Eq. ( I37> can be solved for 
Eb and recovers the universal form 



Eb(B) = - 



m[a(B)Y 



(39) 



of the binding energy. 



Universal bound states 



1. Universal binding energy 

We shall study the binding energy E\, that is determined by 
Eq. j31> in this small region of magnetic field strengths on 
the low field side of the Feshbach resonance. As the binding 
energy vanishes at the resonance position, these studies will 



In accordance with Eq. (I29> . near resonance the highest ex- 
cited bound state is strongly dominated by its component in 
the asymptotic open channel, which is given by 



hir) = 0b^(r) = 



yjlnaiB) r 



(40) 
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at interatomic distances that are large in comparison with 
the van der Waals length. This wave function is ex- 
tended far beyond the outer classical turning point rdasskai - 
a{B) (2Zvdw)^] ^ of the background scattering potential 
Vbgif), with a mean interatomic distance on the order of the 
scattering length; 

<r) = 47r I p-dr \4>h(r)\^ r = a(B)/2. (41) 
Jo 

As a consequence, the coupling between the channels in the 
close vicinity of a Feshbach resonance can be treated, to an 
excellent approximation, as a perturbation of the background 
scattering potential Vbo(r). In fact, to describe the universal 
low energy scattering properties as well as the highest excited 
vibrational bound state, the whole potential matrix can be re- 
placed by a single channel potential V(B, r) that recovers the 
correct magnetic field dependence of the scattering length and 
the length scale associated with the long range van der Waals 
interaction between the atoms. This involves the replacement 
of the two-body Hamiltonian matrix in Eq. (|2jl by the single 
channel Hamiltonian 

H2B = V2 + y(B,r). (42) 

m 

Figure 0] shows a comparison between the component cf)^ of 
the highest excited vibrational bound state as determined from 
a full multichannel calculation, and from a single channel 
Hamiltonian with the background scattering potential modi- 
fied in such a way that it recovers the scattering length of the 
full multichannel Hamiltonian. 

When the magnetic field strength approaches the resonance 
position from below, the spatial extent of (f>h{r) becomes infi- 
nite and the bound state wave function becomes degenerate 
with the zero energy scattering state of two asymptotically 
free atoms in the open channel. Consequently, by sweeping 
the magnetic field strength adiabatically across the resonance 
from negative to positive scattering lengths the zero energy 
scattering state is transferred smoothly into the bound state 
4>h{r). This is the key feature of the adiabatic association of 
molecules in ultracold atomic gases. We note that applying 
the sweep of the magnetic field strength in the direction from 
positive to negative scattering lengths is not suited to asso- 
ciate molecules because there is no energetically accessible 
vibrational bound state on the high field side of the Feshbach 
resonance. Both the association mechanism and its asymme- 
try in the direction of the sweep can be understood purely in 
terms of two-body considerations. 



D. Transition probability 

We shall show in the following how the dynamics of the adi- 
abatic association technique can be described in terms of the 
previous considerations. We assume that the magnetic field 
strength is swept linearly across the resonance position start- 
ing asymptotically far from the resonance on the high field 




10^ lo' lO"* 10^ 

^ BohrJ 

FIG. 4: Radial wave functions associated with the component (j)'^ 
of the highest excited vibrational bound state as determined from a 
full multichannel calculation (solid curve) and from a single channel 
Hamiltonian (dashed curve) with the background scattering potential 
modified in such a way that it recovers the same scattering length 
a{B) = 7200 QBohr as the full multichannel Hamiltonian. The dotted 
curve shows the result of an analogous calculation with the minimal 
two channel Hamiltonian in the appendix. The components (p'^^ of the 
two- and multichannel wave functions agree at interatomic distances 
r that are large in comparison with the van der Waals length of about 
80 flBohr- The minimal two channel description does not account for 
all the bound states of the exact background scattering potential. The 
nodes of the multichannel wave function are therefore not recovered. 
The differences between the single channel wave function and the 
solid and dotted curves result from their different normalisation due 
to the small but still relevant 10 % admixture of the resonance state in 
the multichannel wave functions. We note that the radial coordinate 
is given on a logarithmic scale. 

side. The two-body Hamiltonian in Eq. Q becomes explic- 
itly time dependent through its dependence on the magnetic 
field strength, i.e. H2B = H2B.it)- The probability for the adi- 
abatic association of a pair of atoms in the initial state I*!*,) is 
determined in terms of the two-body time evolution operator 
U2B(tf,ti) by: 

Pfi^\{MBf)\U2B(tf,tiWi}f . (43) 

Here f, and tf are the initial and final times before and after 
the linear ramp of the magnetic field strength, respectively, 
while 0b(fi /) is the highest excited vibrational bound state of 
the two-body Hamiltonian //2B(f/) at the final magnetic field 
strength B(f/). The time evolution operator is determined by 
the Schrodinger equation: 

d 

in-U2B(t, ti) = H2Bm2B(t, f/). (44) 

ot 

The transition amplitude in Eq. J43> for the minimal two chan- 
nel Hamiltonian in the appendix can be obtained with meth- 
ods similar to those applied to the determination of the single 
channel two-body time evolution operator in Ref . Il45tl . A tran- 
sition amplitude similar to that in Eq. ( I43t serves as an input to 
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the microscopic quantum dynamics approach to the associa- 
tion of molecules in a trapped dilute Bose-Einstein condensate 
that is presented in Subsection llll Al 

We shall use Eq. J43> to determine the probability for the 
association of two atoms in the ground state of a large box of 
volume 'V that is later taken to infinity. In the homogeneous 
limit the appropriate initial state in Eq. j43t is given by the 
product state 



l*P/) = |0,bg) 



(45) 



where |0) is the isotropic zero momentum plane wave of the 
relative motion of two atoms in free space and |bg) denotes 
their internal state in the asymptotic open channel. The factor 
-sjilnti)^ in Eq. ( 145 > provides the appropriate normalisa- 
tion. In the homogeneous limit the product P ffV is indepen- 
dent of the volume 'V. 

Figure|5]shows the product PffV as a function of the initial 
time f, for a linear sweep of the magnetic field strength across 
the 100.74 mT Feshbach resonance of ^^Rb with a ramp speed 
of 0.1 mT/ms. The calculations were performed on the basis 
of Eq. ( I43> with the two channel approach of the appendix 
(solid curve) and with a single channel Hamiltonian (dashed 
curve) that properly accounts for the width of the resonance. 
Although the detailed time evolution is slightly different in 
the single and two channel case the transition probabilities are 
virtually equal once the ramp starts sufficiently far outside the 
width of the resonance. We shall show in Subsection lll El how 
this independence of the final molecular population from the 
details of the two-body description can be derived in a system- 
atic way. This derivation will reveal that in the limits f, — > -oo 
and r/ — > oo the product PffV depends only on the atomic 
mass m, the background scattering length Obg, the width of the 
resonance (AB), and the ramp speed. 



E. Landau-Zener approach 

An intuitive extension of the two-body dynamics to the adi- 
abatic association of molecules in a dilute Bose-Einstein con- 
densate has been developed by Mies et al. IsSl . This ap- 
proach is based on the time dependent two-body Schrodinger 
equation 



m-^>{t) = i/2B(f)*P(0 

Ot 



(46) 



in the spherically symmetric harmonic potential of an opti- 
cal atom trap 146]. The confining atom trap modifies the 
potentials in Eq. (|3 to Vhgif) — > Vheif) + 57*^^0'"^ 
Vc\{B,r) — > Va(.B,r) + jj^^^^r^, where Who is the angular 
trap frequency. To extend the binary dynamics to the associa- 
tion of molecules in a dilute Bose-Einstein condensate, Mies 
et al. have formulated Eq. i46t in terms of a basis set expan- 
sion with respect to the single channel energy states in the 
open channel and the closed channel strongly coupled to it. In 
this subsection, we shall develop an improved version of this 
approach. 
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FIG. 5: The product P/i'V as a function of the initial time f, for two 
*^Rb atoms in the ground state of a large box with the volume "V. 
The solid curve shows results based on the two channel approach of 
the appendix, while the dashed curve corresponds to an analogous 
calculation with a single channel Hamiltonian that properly accounts 
for the width of the resonance. The final time in these calculations 
is tf = 80 ytis; tj = corresponds to the initial time at which the 
linear ramp of the magnetic field strength starts at the position of the 
Feshbach resonance. In the case of ?, < the linear ramp crosses 
the resonance position from above. The horizontal arrow indicates 
the asymptotic Landau-Zener prediction for Pfj'V as obtained from 
Eq. <61> in Subsection lll El 



1. Two-body configuration interaction approach 



We label the spherically symmetric energy states associated 
with the background scattering by the vibrational quantum 
numbers v - 0,1,2, These energy states fulfil the sta- 
tionary Schrodinger equation 



m 



i,(r) = Ey4>„{r). 



(47) 



In contrast to the free space continuum energy states (p^p\r) in 
Eq. (|8}, the vibrational trap states 0i.(r) are confined in space 
and we shall assume them to be unit normalised. In the atom 
trap the dissociation threshold energy in the open channel is 
given by E^. In analogy to Subsection lll CI we will label the 
vibrational bound states (p^ir) below this threshold by nega- 
tive quantum numbers v = -1, -2, -3, .... For realistic atom 
traps and background scattering potentials, we can presuppose 
the condition lis-il » fiajho to be fulfilled. As a consequence, 
the vibrational bound states below the dissociation threshold 
are hardly modified by the trapping potential. We shall as- 
sume furthermore that the configuration of the atoms in the 
closed channel that is strongly coupled to the open channel, 
is restricted to the resonance state 0res('")- This assumption is 
analogous to the pole approximation in Eq. ( I15l l. The basis 
set expansions for the open channel and the closed channel 
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components of *P(f) are then given by 

I' 



(48) 



(49) 



respectively. The expansion coefficients Cy(f) and CresiO are 
determined by Eq. i46i in terms of the equivalent infinite set 
of coupled configuration interaction (CI) equations: 

inC,(t) = E,Q.(t) + {'/>AW\<Pr..)Q,At), (50) 

inC,,M = £res(B(f))Cres(f) + J]{'f>r.s\W\(P„)C,it). (51) 

V 

These equations determine the exact dynamics of the adiabatic 
association of two atoms in the confining potential of an atom 
ti-ap. 

2. Many-body configuration interaction approach 

There are two essential phenomena that should be taken 
into account to extend the two-body CI approach to the many- 
body physics of the adiabatic association of molecules in a 
dilute Bose-Einstein condensate. First, the many-body mean 
field interactions can cause the size of the atom cloud to oc- 
cupy a much larger volume than a single atom in the harmonic 
trapping potential. This size is determined by the nonlinear- 
ity parameter k\,^ - Na\,g/lho of the Gross-Pitaevskii equation 
i47il . Here is the number of atoms and Zho = V^/(wjJho) 
is the oscillator length of the atom trap. We shall assume in 
the following that the Thomas Fermi condition k\,g » 1 is 
fulfilled. This directly implies that the extension of the atom 
cloud is characterised by the Thomas Fermi radius: 



'tF = /ho (15 fcbg) 



1/5 



(52) 



This radius determines the mean kinetic energy per atom 



{£'kin) , 5 / /ho 



In 



TF 



1.2683 /h 



(53) 



The physical intuition underlying the Thomas Fermi limit re- 
lies upon the observation that the mean field potential cancels 
the trap potential at distances smaller than the Thomas Fermi 
radius. Under these assumptions the single atoms experience 
an eff'ective flat potential rather than the harmonic potential of 
the atom trap. Motivated by this physical intuition, Mies et al. 
replaced the vibrational trap states in the matrix elements in 
Eqs. (15 0> and (15 U by those that correspond to a spherical box 
with a zero point energy of {E\^a)IN. The radius of this box is 
then determined by: 



/bo 



h 



1/2 



a/'" ( 1.2683 /hJ 



(54) 



The second many-body phenomenon taken into account by 
Mies et al. is the macroscopic occupation of the lowest en- 
ergy mode corresponding to the Bose-Einstein condensate. 



Neglecting the curve crossing between E^es(B) and the energy 
£■_! of the highest excited bound state of the background scat- 
tering potential in Fig. |2l in a downward ramp of the mag- 
netic field strength across the Feshbach resonance the prevail- 
ing asymptotic transition involves only the condensate mode 
and the resonance state. As each atom has N - I atoms to 
interact with, the element of the potential matrix that is asso- 
ciated with this prevailing transition is enhanced by a factor of 

Combining these results leads to a two level Rabi flopping 
model for the adiabatic association of molecules in a Bose- 
Einstein condensate: 



inCo(t) = £oCo(f) + ^SQ*C,-es(f), 
Here the Rabi frequency is given by 



Q. = 2 VA^ - 1 



1 



<0res|W|<^J,^'). 



\ 3 'l 



3 

box 



(55) 
(56) 

(57) 



We note that we have quoted the Rabi frequency in terms of 
the free space zero energy background scattering state <^Q^'(r) 

as introduced in Eq. (|8}. The factor -yjilnh)^/ (^/^o^) P'^^" 
vides the proper normalisation that accounts for the finite vol- 
ume 'V = Y^box °^ '■^^ accordance with Eq. ( I25l l, the 
matrix element (^i-eslW'I'^o^^) in Eq. (j^) can be expressed in 
terms of the resonance width (AB), the slope [^^(^res)] of 
the resonance and the background scattering length dbg- 

When the magnetic field strength is swept linearly across 
the Feshbach resonance the energy of the resonance state 
changes linearly in time: 



£res(B(f)) = Eo + 



dEr^, 



dB 



-(Bi-es) 



-T-(fres) 

dt 



(f-f,-es). (58) 



Here fres is the time at which the energy of the resonance state 
crosses the dissociation threshold energy of the open channel, 
i.e. B(fies) = fires and £'ies(fiies) = £^0- Under the assump- 
tion of a linear Feshbach resonance crossing with the initial 
populations |Co(f,)P = 1 and |Cres(f,)P = 0, the final popula- 
tions |Co(f/)P and |Cres(f/)P can be determined by the Landau- 
Zener formulae: 



|Co(f/)P 

IQes(f/)l' 



1 -g-Z'^'SLZ 



(59) 
(60) 



in the limits f. 



-oo and tf — > oo. Derivation of the asymp- 



totic Landau-Zener populations requires much tedious calcu- 
lation. Given the known general form of the exponent 5\jz, 
however, a short calculation using Eq. M5\ reveals its simple 
dependence on the background scattering length Obg, the res- 
onance width (AB) and the ramp speed |^(?res)|: 



4|^(fires)||f(f.-e,s) 



(A^- l)47r;!|flbg||AB| 

^«t|f(fres)| 



(61) 
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Although Eq. i61\ can be derived on the basis of the simple 
two level Rabi flopping model in Eqs. i55\ and j56> . the two- 
body (A^ = 2) Landau-Zener prediction is accurate even in ap- 
plications to transition probabilities that include a continuum 
of energy levels above the dissociation threshold (see Fig.|5j- 
In fact, the Landau-Zener coefficient for the asymptotic popu- 
lation of the resonance state 0ies('") can be derived rigorously 
for an arbitrary number of linear curve crossings associated 
with a quasi continuum of two-body energy states We 
note, however, that despite the universality of the asymptotic 
populations, a two level model is not suited to provide an ad- 
equate description of the intermediate states and their inter- 
mediate populations. From this viewpoint, the agreement be- 
tween the asymptotic two- and multilevel descriptions is co- 
incidental. 

When applied to a gas of many atoms the linear two level 
Rabi flopping model in Eqs. ( I55t and ( I56t yields analytic pre- 
dictions on the efficiency of molecular production. The treat- 
ment of the Bose enhancement in the linear model, however, 
does not account for the depletion of the condensate mode 
in the course of the adiabatic association. This depletion can 
be accounted for in a straightforward way by replacing the 
initial number of condensate atoms by the actual number 
N\Co{t)\^ in the enhancement factor This inclusion of the de- 
pletion modifies the linear Rabi flopping model in Eqs. i55\ 
and ( I56> to the nonlinear dynamic equations: 

mCoit) = EoCoit) + i^Q*|Co(?)|C,es(f), (62) 

inCr,,(t) = £,e.s(B(0)Cres(f) + ^ftQ|Co(0|Co(f). (63) 

At final times f/- — » oo we shall interpret A^|C()(f/)p as the 
number of atoms in the remnant Bose-Einstein condensate and 
A^|Ci-es(f/)P as the number of atoms converted into molecules. 
The number of diatomic molecules produced in the adiabatic 
association is then given by A^|Cres(f/)|^/2. We show in Sub- 
section IIII CI that Eqs. i62i and J63l l significantly improve 
the predictions of the asymptotic Landau-Zener formulae in 
Eqs. ( 15 9> and (I60> . in particular when the molecular produc- 
tion begins to saturate. 

F. Dissociation of molecules 

We note that the chemical bond of the diatomic molecules 
shifts the atomic spectral lines. This line shift can prevent 
the bound atoms from scattering light of probe lasers even in 
the case of the very weakly bound Feshbach molecules that 
are produced with the adiabatic association technique. As a 
consequence, many present day experimental molecular de- 
tection schemes rely upon the spatial separation of bound and 
free atoms in the cloud and the subsequent dissociation of the 
molecules (cf., e.g., Ref. 1 10]). The highest excited vibrational 
molecular bound state can be dissociated by crossing the Fes- 
hbach resonance from positive to negative scattering lengths 
which corresponds to an upward ramp in the case of the 100 
mT Feshbach resonance of'^^Rb. The crossing of the Fesh- 
bach resonance transfers the bound molecules into correlated 



pairs of atoms which can have a comparatively high relative 
velocity, depending on the ramp speed. In ffie following we 
shall characterise the dissociation energy spectra for the 100 
mT Feshbach resonance of '^^Rb, under the assumption that 
many-body phenomena can be neglected. We expect this as- 
sumption to be accurate because the final continuum states of 
the molecular fragments have a low occupancy. Consequently, 
phenomena related to Bose enhancement should be neghgible. 

1. General determination of dissociation energy spectra 

We shall consider linear upward ramps from the initial mag- 
netic field strength B, at time f, across the Feshbach reso- 
nance to the final field strength Bf at time tf. Starting from 
the highest excited multichannel vibrational molecular bound 
state <pb{Bi) at ffie magnetic field strength B,, the state of a 
pair of atoms at any time t is determined in terms of the time 
evolution operator in Eq. ll44t by: 

^^(O^ f/2B(f,f,)0b(B,). (64) 

The dissociation energy spectrum is usually measured in a 
time of flight experiment that allows the fragments to evolve 
freely after the final time tf of the ramp. At any time t after 
tf the two-body Hamiltonian is stationary, and the time evo- 
lution operator can be factorised into its contribution of the 
linear ramp of the magnetic field strength and a part ffiat de- 
scribes the subsequent relative motion of a pair of atoms: 

f/2B(f, ti) = f/2B(f - f/)f/2B(f/, ti). (65) 

Here U^ait - tf) is determined in terms of ffie stationary two- 
body Hamiltonian H2B.{Bf) at the final magnetic field strength 
By by: 

t/2B(f-fy) = e-'«^B(B/»-//)/fi_ (66) 

We shall insert Eq. \65\ into Eq. ( I64> and represent the time 
evolution operator of the relative motion of a pair of atoms 
after the time tf in terms of the complete set of multichannel 
energy states at the final magnetic field strength Bf. The state 
*l'(f) can ffien be decomposed as 

"^(t) = *l'free(0 + *I'bound(0, (67) 

where *l'free(f) describes the dissociation into two asymptoti- 
cally free fragments, while *l'bound(0 describes those events in 
which the bound state <pb(Bi) is transferred into more tightly 
bound multichannel molecular vibrational states at the final 
magnetic field strength Bf. It is the continuum part *l'free(0 
of the wave function ^'(t) that determines the measurable 
dissociation energy spectrum. In terms of the continuum 
multichannel energy states 0p(B/) at the final magnetic field 
strength B f this is given by 

*Pfree(0 = j dp(f>piBf)e-'^^'-'f^'\(/>p{Bf)\U2B(tf,ti)\<Pb{Bi))- 

(68) 
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Here E - p^/m is the energy of the relative motion of the frag- 
ments that corresponds to their relative momentum p. From 
Eq. (I68> we deduce the probability of detecting a pair of atoms 
with a relative energy between E and E + dE to be: 

n(E)dE^p^dp j t/Qp |<0p(B/)|f/2B(f/,f,)l0b(B,))p. (69) 
Here dSlp denotes the angular component of dp. 

2. Exact dissociation energy spectra 
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FIG. 6: Dissociation spectra of the highest excited vibrational bound 
state of *'^Rb2 as a function of the relative energy of the fragments. 
The speeds of the upward ramps across the 100 mT Feshbach reso- 
nance were varied between 0. 1 mT/ms (dashed dotted curve) and 1 
mT/ms (solid curve). The spectra are rather insensitive to the range 
of magnetic field strengths covered by the ramp; the solid curve in- 
dicates a 1 mT/ms ramp with 6, = 100.7 mT and Bf = 100.78 mT, 
while the data indicated by the dots on top of the solid curve cor- 
respond to a calculation for a 1 mT/ms ramp with initial and final 
magnetic field strengths that are half way closer to the resonance po- 
sition of Bo = 100.74 mT. We note that the energies are given on a 
logarithmic scale. 

Figure|6]shows the spectral density n{E) for different ramp 
speeds as obtained from an exact solution of the Schrodinger 
equation i44l with the low energy two channel Hamiltonian in 
the appendix. Although the low energy two-body Hamiltonian 
supports a comparatively tightly bound state at the high field 
side of the Feshbach resonance (see Fig. O, the calculations 
do not indicate any transfer into this state for the reaUstic ramp 
speeds under consideration. Within this range of ramp speeds 
between 0.1 and 1 mT/ms, the spectra cover a broad range 
of energies of the relative motion of a pair of atoms that is 
much larger than the typical energy spread of a Bose-Einstein 
condensate. Consequently, the molecular fragments will be 
detected as a burst of correlated pairs of atoms moving radially 
outward from the original position of the molecular cloud. 



In the following, we shall study in more detail the depen- 
dencies of the dissociation spectra in Fig.|6lon the five physi- 
cal parameters of a Feshbach resonance (cf. Subsection lU Bl i. 
which completely characterise the resonance enhanced low 
energy collision physics. To this end, we shall neglect the 
vibrational bound states of the background scattering poten- 
tial and reformulate Eq. ( I69I I in terms of the energy states in 
the absence of interchannel coupling, in the formal asymp- 
totic limits f; — > -oo and tf — > oo. The idealising assumption 
that the background scattering potential does not support any 
vibrational bound states implies that the two channel Hamilto- 
nian in Eq. Q supports just the initial vibrational bound state 
<ph{Bi) of Eq. \69\ . In the formal limit f, — > -oo, i.e. when 
the linear ramp of the magnetic field strength starts asymptot- 
ically far from the resonance position, (pbiBj) is then identical 
to the closed channel resonance state (p^es- Furthermore, in 
accordance with Eqs. ( I16l l. ( I17> and ( I19t . the final continuum 
state (f>p(B f) of Eq. i69\ is transferred into the energy state 0p^* 
of the background scattering [cf. Eq. (|8}], in the formal limit 
tf — > oo. Under these assumptions, Eq. ( I69l l can be reformu- 
lated to be 



n(E)dE 



P'dp 



bg|t/2B(f/,f,)l0re,s,cl>|'. (70) 



Here |bg) and |cl) denote the internal states of an atom pair in 
the asymptotic open channel and the closed channel strongly 
coupled to it, respectively. The asymptotic dissociation spec- 
trum in Eq. ilOl can then be determined from the results of 
Ref. i49ll in analogy to the derivation of the Landau-Zener 
formulae of Subsection lU El This yields the analytic result for 
the spectral density [501 : 



oE 




(71) 



We note that the asymptotic energy density of the disso- 
ciation spectrum in Eq. illl depends, like the Landau-Zener 
coefficient in Eq. (16 U . only on those physical parameters that 
determine the universal low energy scattering properties in the 
close vicinity of a Feshbach resonance (cf. Subsection 111 C> . 
while the slope [^^(^res)] of the resonance and the van der 
Waals dispersion coefficient C(, do not contribute to Eq. J71l i. 
This corroborates the observation in Fig. [S] that the relevant 
physics occurs in a small region of magnetic field strengths 
in which the modulus of the scattering length by far exceeds 
all the other length scales set by the binary interactions. In 
our applications to the 100 mT Feshbach resonance of ^^Rb, 
the results obtained from Eq. i7l\ are virtually indistinguish- 
able from those of the exact calculations in Fig.|6l We note, 
however, that Eq. (I71> may not be applicable when the ramp 
of the magnetic field strength starts in the close vicinity of a 
Feshbach resonance. This situation may occur, for instance, 
in experimental applications involv ing the broad [(AB) =1.1 
mT] resonance of '^^Rb at 15.5 mT jsj. 
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4. Mean kinetic energy 



III. ADIABATIC ASSOCIATION OF MOLECULES IN A 
TRAPPED BOSE-EINSTEIN CONDENSATE 



In the following, we shall study the mean kinetic energies 



(^kin) 



2X 



dE En(E) 



(72) 



of the molecular fragments after the dissociation. These ener- 
gies characterise the speed of expansion of the gas of molec- 
ular fragments before the detection in related experiments 
lIToll . We note that the kinetic energy of a single atom is 
E\dn - P^/C^m), which is half the energy of the relative motion 
of a pair Equations J71> and (I72> allow us to represent (/Skin) 
in terms of physical parameters of a Feshbach resonance, of 
the ramp speed and of Euler's F function: 



(^kin) - ^ 



S2 /!|f(fn 



|AB| 



2/3 



F(2/3). (73) 



We expect this prediction to be accurate provided that the ini- 
tial and final magnetic field strengths of the linear ramp are 
asymptotically far from the resonance position. 




0.4 0.6 
Ramp speed [mT/ms] 

FIG. 7: Mean single particle kinetic energies of the molecular frag- 
ments associated with the exact dissociation spectra in Fig.|6|(circles) 
in dependence on the ramp speed. The solid curve has been obtained 
from the asymptotic prediction in Eq. J73> . The single particle kinetic 
energy is half the energy of the relative motion of the fragments in 

Fig.m 

Figure shows the strong dependence of the mean single 
particle kinetic energies on the ramp speed. The circles indi- 
cate the mean energies of the exact dissociation energy spec- 
tra in Fig. |6l while the solid curve has been obtained from 
Eq. i73l . Both approaches lead to virtually the same pre- 
dictions for the realistic ramp speeds and physical parame- 
ters of the 100 mT Feshbach resonance of ^^Rb. Figure 
also reveals that ramp speeds of less than about 0.03 mT/ms 
would be required to suppress the mean single particle ener- 
gies of the fragments to a range below (£'kin)/^B = 100 nK 
(fefi = 1.3806503 X 10-23 j/K), which would be close to the 
typical energy spread in a Bose-Einstein condensate |47i|. 



The results in Subsection III Dl allow us to rigorously treat 
the two-body physics of the adiabatic association. In this sec- 
tion we shall study the many-body physics of the molecular 
production in a trapped dilute Bose-Einstein condensate. In- 
homogeneous Bose-Einstein condensates are subject to a rich 
spectrum of collective energy modes that depend sensitively 
on the number of condensate atoms and their binary interac- 
tions, as well as on the confining atom trap. The crossing 
of a singularity of the scattering length leads to a dramatic 
change of the intermediate collision physics that transfers a 
substantial fraction of the initially coherent atomic cloud into 
strongly correlated pairs of atoms in the highest excited vibra- 
tional bound state. This violent dynamics can be expected to 
couple all energy modes. A proper description of the com- 
plex interplay between the macroscopic collective behaviour 
and the microscopic binary collision physics requires a full 
many-body treatment. We shall provide an appropriate many- 
body description of the adiabatic association of molecules in 
a trapped dilute Bose-Einstein condensate and compare it to 
previous theoretical predictions at different levels of approxi- 
mation. This will highlight in what parameter regimes approx- 
imate descriptions are valid and to which physical observables 
they are applicable, as well as when we expect them to break 
down. 



A. Microscopic quantum dynamics approach 

The general many-body approach to the dynamics of 
atomic gases that we shall apply has been derived in Refs. 1 2^ 
|4^ . This approach has been applied previously to several 
different physical situations that involve the production of 
correlated pairs of atoms in four wave mixing experiments 
with Bose-Einstein condensates |26], the determination of the 
mean field energy associated with three-body collisions f52|], 
and the dynamics of atom molecule coherence 1 45, 53] as well 
as Feshbach resonance crossing experiments with degener- 
ate Bose gases of ^^Rb and ^^Na atoms l54ll . The underly- 
ing method of cumulants (25 Issll is based on the exact time 
dependent many-body Schrodinger equation. As the general 
technique has been derived before, we shall outline only those 
details specific to the adiabatic association of molecules in a 
Bose-Einstein condensate. 



1. Multichannel approach 

We shall formulate the approach in terms of the full multi- 
channel many-body Hamiltonian, which couples all the inter- 
nal hyperfine states of the single atoms. We label the quan- 
tum numbers associated with the single atom energy states by 
Greek indices. In our applications to the low energy scatter- 
ing in the vicinity of the 100 mT Feshbach resonance of ^^Rb, 
these indices are sufficiently characterised by the total angular 
momentum quantum number F and its orientation quantum 
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number nip. In its second quantised form, the full many-body Hamiltonian then reads: 
I 



" ^ fJVKA ^ 



(74) 



r 



Here the single particle annihilation and creation field opera- 
tors satisfy the Bose commutation rules: 



[ifffjix), iff'liy)] = 6^y6(x - y), 
[tA;.(x),tAv(y)] = 0. 



(75) 
(76) 



Furthermore, Hj^iB) is the one body Hamiltonian associated 
with the internal atomic state yU that contains the kinetic energy 
of the atom, the external potential of the optical atom trap, and 
the internal hyperfine energy E'^(E): 



2m 



+ V, 



trap 



+ <(B). 



(77) 



The hyperfine energy depends on the magnetic field strength 
through the Zeeman effect. The microscopic binary potential 
^l/iv),|«^l(r) in Eq. is associated with the asymptotic in- 
coming and outgoing binary scattering channels that can be 
labelled by the pairs of internal atomic quantum numbers {kA} 
and {yuv), respectively. We note that in this general formula- 
tion of the many-body Hamiltonian all potentials associated 
with the asymptotic binary scattering channels are chosen to 
vanish at infinite interatomic distances. 

All physical properties of a gas of atoms are determined 
by correlation functions, i.e. quantum expectation values (de- 
noted by (■),) of normal ordered products of field operators 
for the quantum state at time t. The correlation functions 



J 



of main interest in the adiabatic association of molecules in- 
volve the atomic mean field ((/(^(x)),, the anomalous average 
(i/'y(y)i/'^(x)),, and the one-body density matrix (i^i(y)i/'^(x)),. 
The dynamics of the correlation functions is determined by 
the many-body Schrodinger equation through an infinite hier- 
archy of coupled dynamic equations. The cumulant approach 
consists in transforming the coupled set of dynamic equations 
for the correlation functions into an equivalent but more prac- 
tical infinite set of dynamic equations for noncommutative cu- 
mulants |26, 45]. The transformed equations of motion for the 
cumulants can be truncated, at any desired degree of accuracy, 
in accordance with Wick's theorem of statistical mechanics 
js^l- The noncommutative cumulants that we shall consider 
in the following are the atomic mean field 



'¥^{x,t)^{Mx)}„ 

the pair function 

%yix, y, t) = {ilf,(y)ilf^(x)), - >P^(x, f)*I'v(y, 0, 
and the density matrix of the noncondensed fraction 



(78) 



(79) 



(80) 



r;.v(x, y, t) = {^l{y)ilf,{x)), - *I';.(x^)r(y, f). 
In the first order cumulant approach 12a |451 that we shall 
apply in this paper the atomic mean field and the pair function 
are determined by the coupled dynamic equations: 



i^j^.i^, t) ^H]^{B)^,{x, t) + J]fdy X(y, OV,;.v.),|.,i(x - y) [(D,,(x, y, t) + ^^.(x, ^^^^(y, 0] (81) 

vkA 

d 

in-%Ax, y, = X ["ZlWA](B)^'<^iX' y- 1) + ^U-v)>i)(x - tm(j, t)] . (82) 

I 



kA 



Here H2b{B) is the Hamiltonian matrix of two atoms whose 
matrix elements 

<W>^)('S) = [-f^f (B) + Hl^'iBm^dAy + y,^v),k^l (83) 

involve all incoming and outgoing asymptotic binary scatter- 
ing channels associated with the pairs of atomic indices {kA} 
and (yuv), respectively. Given the solution of Eqs. j81> and 
(I82> . the density matrix of the noncondensed fraction is deter- 



mined in terms of the pair function by ll45ll 

r^v(x, x', f) = 2 J fl'y y, miJx', y, t). (84) 



The general first order dynamic equations (I81> . (I82> and (|84| 
not only strictly conserve the total number of atoms 



Yj [dx [r^:^(x,x,t) + \^'^(x,t)f■] (85) 
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at all times, but the explicit form of Eq. (I84> also ensures the 
crucial property of positivity of the one-body density matrix. 

Equation ( I84> reveals that the density of the noncondensed 
fraction stems directly from the pair function. Consequently, 
the build up of pair correlations is the main source of atom 
loss from a Bose-Einstein condensate. The energy delivered 
by a time dependent homogeneous magnetic field can trans- 
fer a pair of condensate atoms either into the highest excited 
vibrational molecular bound state or into the quasi continuum 
of two-body energy states above the dissociation threshold. 
Not only the molecular association, but also the production of 
correlated pairs of atoms in the scattering continuum has been 
observed recently 12J as a burst of atoms ejected from the rem- 
nant condensate. As a time dependent homogeneous magnetic 



field delivers energy but no momentum to the gas, the centres 
of mass of all correlated bound and free pairs of atoms have 
the same momentum distribution as the initial Bose-Einstein 
condensate. In this sense, the molecules produced in the adia- 
batic association may be considered as a degenerate quantum 
gas. The number A^b of diatomic molecules in the state 0b is 
determined by counting the overlap of each pair of atoms with 
the multichannel molecular bound state |57]. This relates the 
number of molecules to the two-body correlation function 



Gfl/x, y; x', y') = <^l(x')^l(y')<Av(y)«A^(x)> (86) 



by 
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A^b 
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(87) 



r 



Here the spatial integration variables can be interpreted in 
terms of two-body centre of mass and relative coordinates 
R - (x + y)/2 and r = x - y, respectively. The number of 
correlated pairs of atoms in the scattering continuum can be 
deduced from the two-body correlation function in a similar 
way 1 45] with the multichannel bound state wave function re- 
placed by the continuum states. We note that Eq. (|^) neither 



assumes any particular class of many-body states nor any ap- 
proximation to the many-body Schrodinger equation. By ex- 
panding the two-body correlation function in Eq. j87t into cu- 
mulants and truncating the expansion in accordance with the 
first order cumulant approach, the density of bound molecules 
can be represented in terms of a molecular mean field 



J 



^b(R, = ^ f^*" K''"'^'")]* [•^'^''(^' r, + (r + ^, f) 1'v (r - ^, f) 



(88) 



r 



Here we have introduced the centre of mass and relative co- 
ordinates R and r and represented the pair function in terms 
of these variables. The molecular mean field determines the 
density of diatomic molecules in the state 0b by «b(R> - 
|*l'b(R, 01^- The molecular mean field as well as the fraction 
of pairs of correlated atoms in the scattering continuum are de- 
termined completely by the solution of the coupled equations 
n\ and (|83. 



2. Two channel approach 

The general form of the two-body Hamiltonian in Eq. dSSl 
with a realistic potential matrix allows us to describe the bi 
nary collision physics over a wide range of energies and mag 
netic field strengths, as indicated in Fig.[0 As the present ap 



plications involve only the adiabatic association of molecules 
in the vicinity of the 100 mT Feshbach resonance of **^Rb, we 
can restrict the asymptotic binary scattering channels to those 
that we have identified in Section [II] We shall thus insert the 
two channel description of Section Hll into Eqs. ( I81> and i&2\ 
and perform the pole approximation of Subsection lllBI The 
relevant potentials then involve the background scattering po- 
tential Vh^{r) and the off diagonal matrix element W{r) be- 
tween the open channel and the closed channel strongly cou- 
pled to it. In accordance with the pole approximation, the 
only configuration of the atoms in this closed channel is the di- 
atomic resonance state 0res('")- Consequently, the atomic mean 
field is restricted to its component in the {F - l,mf - +1) 
state, which we shall denote simply by *l'(x, t). According to 
Eq. ( 18 1> , the associated mean field dynamic equation is then 
given by 
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"Vix, t) + Jdy >P*(y, f) (W(|x - y|)<Dei(x, y, t) + V^,(\x - y|) [Obg(x, y, t) + >y(x, tmy, t)]) 



(89) 



r 



The pair function has a component in the open channel and in Eq. J82l i. their coupled dynamic equations read: 
the closed channel strongly coupled to it. In accordance with 



d_l a)bg(x,y,f) 
dt\ 1>ci(x,y,f) 



rriB 3>bg(x,y,f) 
^trap^r^l Od(x,y,f) 
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¥(x,f)»P(y,f)- 



(90) 



Here H^^t) is the general two channel two-body Hamiltonian 
[cf. Eq. that includes the centre of mass kinetic energy of 
a pair of atoms as well as the confining harmonic potential of 
the atom trap. In the two channel formulation of the first order 
microscopic quantum dynamics approach, Eqs. j89> and j90> 
completely determine the dynamics of the coherent atomic 
condensate and the fraction of correlated pairs of atoms. 

For the purpose of numerical convenience, we shall solve 
the inhomogeneous linear Schrodinger equation ( I90> formally 
in terms of the complete two-body time evolution operator 



t/y.2p(f, r) [cf. Eq. J44H that accounts for the centre of mass 
motion of the atoms and for the confining trap potential. We 
shall then insert the solution into the mean field equation j89> 
to eliminate the pair function. We shall assume that at the ini- 
tial time f, at the start of the ramp the gas is well described 
by a dilute zero temperature Bose-Einstein condensate, so 
that initial two-body correlations are negligible. The result- 
ing dynamic equation for the atomic mean field can then be 
expressed in terms of the time dependent two-body transition 
matrix 



V(t)S{t - r) + -e{t - T)V{t)uZ^{t, t)V(t) 



|bg). 



(91) 
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which involves the potential matrix 



V(f) = 



w 



w \ 



(92) 



and the two-body time evolution operator U'^^^{t, t). Here |bg) 
denotes the internal state of a pair of atoms in the asymptotic 
open scattering channel, and 6{t-T) is the step function, which 



J 



yields unity at f > t and zero elsewhere. A short calculation 
applying the methods derived in Refs. l26l |45l] then shows 
that the elimination of the pair function in Eq. j89t leads to a 
closed dynamic equation for the atomic mean field. Expressed 
in terms of the two-body transition matrix in Eq. J9U . this is 
given by: 



ih-'i'{x,t) 
at 



+ Vu-ap(x) 



^^(X, f) ■ 



dT J dydx'dy' ^\y,t)TZ^ix,yJ-x' ,y ,Tmx' ,Tmy\T). (93) 
I 



In the case of a harmonic trap potential the two-body transi- 
tion matrix in Eq. (I91> factorises into a centre of mass part and 
a contribution from the relative motion of an atom pair In the 
following, we shall assume the confinement of the atom trap 



to be sufficiently weak and the ramp speeds to be sufficiently 
high that the time spent within the width of the resonance is 
much smaller than the trap periods in all spatial directions. 
The centre of mass motion of a pak of atoms then becomes 
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negligible on this time scale. The explicit form of Eq. ( I91> 
also reveals that the time evolution operator of the relative mo- 
tion of two atoms is evaluated only within the spatial range of 
the binary interaction potentials, where it is hardly modified 
by the presence of the atom trap. Consequently, the trap po- 
tential can also be neglected in the two-body time evolution 
operator of the relative motion of an atom pair. Furthermore, 
as the atomic mean field is slowly varying on the length scales 
set by the binary interaction potentials, Eq. ( I93> becomes local 
in the spatial variables and acquires the form |26|,(4JJ: 



at 



-— v" + y,„p(x) 

2m 



"Vix, t) 



->P*(x,f) I c/tT2(x,t)— /i(f,T). (94) 

Jt, OT 

Here we have introduced the coupling function 

h(t, t) = 6(t - T){2nn)\0, bg|y(f)f/2B(f, t)|0, bg>, (95) 

which depends on the time evolution operator U2R{t, t) asso- 
ciated with the relative motion of an atom pair in free space, 
while |0, bg) denotes the zero momentum plane wave state of 
the relative motion of two atoms in the asymptotic open scat- 
tering channel. 



3. Atomic condensate fraction 
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interest in the adiabatic association of molecules is largely in- 
dependent of the details of the implementation of the two- 
body Hamiltonian. The main requirement is that the two- 
body Hamiltonian properly accounts for the magnetic field 
dependence of the scattering length and the near resonant 
binding energy of the highest excited vibrational bound state. 
The differences between a single and a two channel treatment 
(cf. Subsection lll Dl i in the dynamics of the atomic mean field 
in Eq. (|93 are marginal (cf., also, Ref. |54]). The follow- 
ing calculations have, therefore, been performed with a single 
channel two-body Hamiltonian as introduced in Ref. l45ll and 
described in the appendix. In the course of our studies, we 
have solved Eq. (|94} for a variety of linear ramps of the mag- 
netic field strength across the 100 mT Feshbach resonance of 
**^Rb in spherically as well as cylindrically symmetric har- 
monic atom traps. Figure [S] shows the time evolution of the 
number of condensate atoms A^c(0 - J d'n |*l'(x, f)P for linear 
downward ramps with a ramp speed of 0. 1 mT/ms in spher- 
ically symmetric atom traps with the frequencies Vho - 100 
Hz and Vho = 10 Hz. The atomic mean field at the initial time 
f, - has been chosen in each calculation as the ground state 
of the Gross-Pitaevskii equation that corresponds to a dilute 
zero temperature Bose-Einstein condensate with = 50000 
atoms and a scattering length of Obg - 100 OBohi - Figure |8l 
reveals that the loss of condensate atoms mainly occurs dur- 
ing the passage across the Feshbach resonance. In accordance 
with the higher local density, the losses of condensate atoms 
are more pronounced in the comparatively strongly confining 
100 Hz atom trap. The pronounced oscillations immediately 
after the passage across the resonance (f > 500 yus) indicate a 
rapid exchange between the condensed and the noncondensed 
phases of the gas. Exchange on these short time scales sug- 
gests that the crossing of a Feshbach resonance drives the gas 
into a highly excited nonequilibrium state. 

Figure |9] shows snapshots of the condensate density 
nc(x, f) = |*l'(x, f)p as a function of the radial coordinate r - |x| 
for the ramps and trap parameters in Fig.|S]at the final time 
tf = 580 yus. The pronounced spatial variations of the rem- 
nant condensate density indicate the simultaneous occupation 
of many collective energy modes. We note that these highly 
occupied excited modes are described by a coherent classical 
mean field and are therefore distinguished from the pairs of 
correlated atoms in the initially unoccupied scattering contin- 
uum. In accordance with Eq. j84> . correlated pairs of atoms, 
as described by the pair function <l>(f), constitute the density 
of the noncondensed fraction. 



FIG. 8: The time evolution of the condensate fraction (A'c(0 = 
j dx |4'(x, OP) for a linear downward ramp of the magnetic field 
strength across the 100 mT Feshbach resonance of *''Rb. The solid 
curve shows a calculation for a weakly confining spherical atom trap 
with the frequency Vho = 10 Hz, while the dashed curve corresponds 
to a comparatively strongly confining (vho = 100 Hz) spherical atom 
trap. The ramp speed in both calculations is 0. 1 mT/ms and the initial 
number of condensate atoms is N = 50000. The time at which the 
magnetic field strength crosses the resonance position ist = 500 yus. 

As we have shown in Subsections 111 CI 111 Dl and 111 El the 

two-body time evolution of the main physical observables of 



4. Noncondensed fraction 

We shall study in the following how the final noncondensed 
fraction of the gas is distributed among the bound and free en- 
ergy states of an atom pair We shall focus on linear downward 
ramps of the magnetic field strength across the 100 mT Fes- 
hbach resonance of '^^Rb. Given the solution of the dynamic 
equation j94l i for the atomic mean field, the pair function 0(f) 
is determined completely by the inhomogeneous linear two- 
body Schrodinger equation ( I90l l. We shall assume that the 
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FIG. 9: Snapshots of the remnant condensate density after the cross- 
ing of the 100 mT Feshbach resonance of '^'Rb in spherically sym- 
metric atom traps with the frequencies Vho = 100 Hz (a) and Vho = 10 
Hz (b). The densities correspond to the calculations in Fig.|8|at the 
final time of the ramps. 



trap is switched off immediately after the ramp so that the gas 
expands freely in space. Inserting the complete set of two- 
body energy states (cf. Subsection lllBl i at the final magnetic 
field strength Bf into Eq. J84> . and integrating with respect to 
the centre of mass coordinate R - (x + y)/2 yields the number 
of noncondensate atoms |45]: 

A^ne ^JdRjdp \(R,<f>p(Bf)mtf)}f 

+ JdR \(R,MBf)mtf)}\\ (96) 

Here we have presupposed that only the energetically ac- 
cessible highest excited vibrational multichannel bound state 
4>h{Bj) on the low field side of the Feshbach resonance 
(cf. Subsection 111 C> will be populated in a downward ramp. 
We shall assume furthermore that the final magnetic field 
strength is sufficiently far from the Feshbach resonance that 
the gas is weakly interacting; i.e. npeak[fl(fi/)]"' ^ 1, where 
"peak is the peak density of the gas, and a(Bf) is the final 
scattering length. The second, factorised contribution to the 



molecular mean field on the right hand side of Eq. ( I88> can 
then be neglected |45]. Consequently, the bound state con- 
tribution to the noncondensed density on the right hand side 
of Eq. \96\ determines the number of atoms associated to 
molecules. The contribution of the continuum part of the 
two-body energy spectrum yields the burst fraction, composed 
of correlated atoms with a comparatively high relative mo- 
mentum Ipl in initially unoccupied modes |45]. According to 
Eqs. ( I85> and ( I96l l. the number of condensate atoms, atoms 
associated to molecules and burst atoms then add up to the 
total number of atoms. We note that Eqs. ( I85> and ( I96> are 
also applicable in the small region of magnetic field strengths 
in the close vicinity of the Feshbach resonance, in which the 
gas becomes strongly interacting. The dilute gas parameter 
npeak[fl(fi)]^ is then on the order of unity or larger. Under these 
conditions, however, a separation of the gas into bound and 
free atoms is physically meaningless 
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FIG. 10: Density of *'Rb2 molecules in the highest excited vibra- 
tional bound state after the crossing of the 100 mT Feshbach res- 
onance in spherically symmetric atoms traps with the frequencies 
Vho = 100 Hz (a) and Vho = 10 Hz (b). The molecular densities 
correspond to the calculations in Fig.|S|at the final time of the ramps. 



Figure ^| shows the densities of **^Rb2 molecules in the 
highest excited vibrational molecular bound state at the fi- 
nal magnetic field strength of the linear downward ramps 
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across the 100 mT Feshbach resonance, which are described 
in Fig. |S] The densities have been obtained from the molec- 
ular mean field *Fb(R, ?/) in Eq. (|88j through nb(R, f/) - 
\^h(R, tf)\^ with the bound state wave function (pbiBf) at the 
final magnetic field strength Bj of the ramps. The spatial ex- 
tent of the molecular clouds roughly corresponds to the size 
of the remnant condensate densities in Fig.|9l In accordance 
with the higher local densities, the density of molecules in the 
tight 100 Hz atom trap is higher than in the 10 Hz trap. 
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FIG. 1 1 : The relative loss of condensate atoms [N^Ui)- Nc(tf )]/Nc(ti) 
in linear upward ramps (circles) as well as downward ramps 
(squares) of the magnetic field strength across the 100 mT Feshbach 
resonance of ^'Rb as a function of the inverse ramp speed. The up- 
ward ramps transfer the condensate atoms into correlated pairs of 
atoms in initially unoccupied excited states, while the downward 
ramps adiabatically associate a substantial fraction of the conden- 
sate to bound molecules in the highest excited vibrational state. The 
diamonds indicate the relative number 2Nf,(tf )/Nc(tj) of atoms asso- 
ciated to molecules in the downward ramps. The calculations were 
performed for a cylindrical atom trap with the axial (radial) frequen- 
cies of Vaxiai = 116 Hz (v„diai = 75 Hz) and N = 140000 atoms. The 
small deviations in the calculated loss data of the remnant atomic 
condensate (filled circles) from an entirely smooth curve are due to 
rapid oscillations in the atomic condensate fraction after the passage 
of the Feshbach resonance (see Fig.|H), which have not entirely de- 
cayed at the final times that we have chosen for the different inverse 
ramp speeds. 

Figure shows systematic studies of the dependence of 
the molecular production, after downward ramps of the mag- 
netic field strength across the 100 mT Feshbach resonance of 
^^Rb, on the inverse ramp speed. The calculations have been 
performed on the basis of the dynamic equation for the 
atomic mean field. In these calculations we have chosen the 
potential of a comparatively tight cylindrically symmetric op- 
tical atom trap with axial and radial frequencies (see Fig.llH 
that resemble those in cuiTent experiments on the production 
of ultracold **^Rb2 molecules in Oxford |59]. The initial state 
has been chosen as the ground state of the Gross-Pitaevskii 
equation with Ndti) - 140000 atoms and a scattering length 
of flbg = 100 flBohr- The squares in Fig. II II indicate the rel- 
ative loss of condensate atoms [Ndtj) - Nc{tf)]/Nciti) in the 



downward ramps, where A^c(0 = f dx \^'(x,t)\^ is the num- 
ber of condensate atoms at time t. As expected from the 
two-body considerations in Subsection 111 El the loss curve 
is a monotonic function of the inverse ramp speed and sat- 
urates when the remnant condensate is completely depleted. 
The number of atoms associated to diatomic molecules in 
the highest excited vibrational bound state 4>h{Bf) is given by 
2Nb(tf) = 2jdR |>l'b(R, where >l'b(R, f/) is the molec- 
ular mean field in Eq. j88> . The dependence of the relative 
molecular association efficiency 2Nbitf)/Nc{ti) on the inverse 
ramp speed in the downward ramps is indicated by the dia- 
monds in Fig. [^and follows the same monotonic trend as 
the loss curve of condensate atoms in the downward ramps 
(squares). The quantitative agreement between both curves 
shows that the missing fraction of condensate atoms is trans- 
ferred into diatomic molecules in the highest excited multi- 
channel vibrational bound state. The excitation of atom pairs 
in continuum energy modes is suppressed. The rapid loss of 
condensate atoms, however, leads to an overall heating of the 
atomic cloud due to the excitation of collective energy modes 
(see Fig.|5}. 

The circles in Fig.llllindicate predictions of the total loss 
of condensate atoms in upward ramps of the magnetic field 
strength across the 100 mT Feshbach resonance of ^^Rb. The 
initial state of the gas as well as the confining atom trap in 
these calculations are the same as in the downward ramps. 
The final magnetic field strength B /, after the passage across 
the resonance in an upward ramp, is on the high field side of 
the resonance position. Figure |5]reveals that there is no ener- 
getically accessible bound state at the magnetic field strength 
B f. Consequently, atoms lost from the condensate are entirely 
transferred into the burst fraction. Figure shows that de- 
spite an overall agreement in the monotonic behaviour of the 
condensate loss with increasing inverse ramp speeds, there is 
a small difference in the saturation region between the curves 
for the different ramp directions. This small deviation be- 
tween the curves would be absent in the two-body descrip- 
tion of Section in] and may indicate phenomena related to the 
coherent nature of the gas. 

B. Two level mean field approach 

One of the most common approaches to the association of 
molecules in dilute Bose-Einstein condensates is based on 
a model many-body Hamiltonian that describes atoms and 
molecules in terms of separate quantum fields. The coupling 
between these quantum fields leads to an exchange between 
the different species. This exchange then serves as a model 
for the molecular production. 



1. Markov approximation to the first order microscopic quantum 
dynamics approach 

In order to reveal the physical meaning of the molecu- 
lar quantum field and its relation to the microscopic quan- 
tum dynamics approach of Subsection IIIIAI we shall derive 
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the mean field equations associated with the two component 
model Hamiltonian on the basis of the first order coupled dy- 
namic equations ( I89> and i90\ that determine the atomic mean 
field and the pair function. We shall represent the pair function 
in terms of the centre of mass and relative coordinates R and 
r, respectively. In accordance with the pole approximation of 
Subsection IllBI the closed channel part of the pair function 
can be factorised into the resonance state 0res('") and a centre 
of mass wave function: 

Oel(R, r, t) = V2,^,es(r)^res(R, 0- (97) 

In the following derivation we shall show that the centre of 
mass wave function *Pres(R, is identical to the mean field 



associated with the molecular quantum field in the two com- 
ponent model Hamiltonian. To this end we shall proceed in a 
way similar to the derivation of the dynamic equation ( I94t for 
the atomic mean field in Subsection llll Al except that we shall 
eliminate only the component <l>bc;(R, r, t) of the pair function 
in the open channel from Eqs. ( I89> and (I90> . The formal so- 
lution of the dynamic equation for <l>bg(R, r, f) [i.e. the first 
component of Eq. j90H can be expressed in terms of the two- 
body time evolution operator f/trap(f - t) that corresponds to 
the stationary diagonal element (bg|//2B|bg) of the full two- 
body Hamiltonian matrix i/2B(0 [cf. Eq. J90H in the open 
channel: 



(Dbg(R,r,f) 



r|f/^s (f-f,)|R',r')<l)bg(R',r',f,) 



: J dR'dr' (R, 

+ ^j^ dT j dR'dr' <R, r\ull^(t - T)Vbg\R', r'mR' + r' 12, t)>I'(R' - r' 12, t) 
+ ^j dTjdR' <R, r\ul%(t - t)W\R', ct>,,,) V2>I',es(R', r). 



(98) 



We note that t/trap(^ ~ ''') includes the centre of mass and the 
relative motion of two atoms as well as the confining potential 
of the atom trap. 

To obtain the mean field equations for ^(x, f) and the cen- 
tre of mass wave function *I're,s(R, in Eq. J97I I. we shall in- 
sert Eq. j98> into Eqs. j89t and j90t and perform the Markov 
approximation. The Markov approximation relies upon the 
assumption that the main contribution to the time integrals 
in Eq. j98> stems from a small region near r = f, in which 
the variation of the functions *I'(R' + r' 12, t), ^{R' - r' 12, t) 
and *I'ies(R', r) is negligible. Under this assumption, the func- 
tions can be evaluated at r = f. A similar argument applies 
to the spatial variation of *I'(R' + r' 12, t), *I'(R' - r' /2, t) and 
*I',e,s(R', t) and leads to the replacements R' — > R and r' ^ 
in these functions. The subsequent formal procedure to derive 
the mean field equations is analogous to the derivation of the 
Gross-Pitaevskii equation in Ref. 1 26] and involves neglecting 
the initial pair function Obg(R', r', f,) and the centre of mass 
motion as well as the confining potential of the atom trap in 



the time evolution operator U^^lpit 
then yields the coupled equations 



t). a short calculation 



in-'i'ix,!) 
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h ^ m ^ 



*P(x, 



+gbgi*i'(x,f)r^(x,f) 

+ gres*i'*(x,0V2»I'res(x,0 



(99) 



for the atomic mean field, and 



/»-'Pres(R,f) = 
Of 



dB 

1 



■ (fires) 



^res(R,0 
[B(f)-Bo]^res(R,f) 

(100) 



+ —g,,,^\R,t) 
V2 

for the molecular mean field S. Here we have assumed the 
optical atom trap to be spherically symmetric. The coupling 
constants gbg and gres can be obtained by performing the spa- 
tial as well as the time integrations in Eq. ( I98> in the limit 
f - f, ^ oo. This yields (cf., also, Ref. j^3l)- 



-dbg 

s,3/2 



(101) 
(102) 



Here abg is the background scattering length and 0Q^\r) is the 
zero energy wave function associated with the background 
scattering (cf. Eq. (|S}). We note that the off diagonal cou- 
phng constant ^res is defined only up to a global phase and 
we shall, therefore, choose it to be real. In accordance with 
Eq. \25\ , gies is then determined by the resonance width (AB), 
the background scattering length abg and the slope of the res- 
onance [^^^ (fires)]- The mean field equations j99t and jlOOt 
as well as the coupling constants in Eqs. dlOU and ( I102> are 
analogous to those that have been applied, for instance, in 
Ref. L6 1.1 to describe the atom loss in ramps of the magnetic 
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field strength across Feshbach resonances in ^^Na. In the ap- 
plications of the mean field equations i99\ and jlOO> . the den- 
sities nc(x, f) = |*l'(x, f)P and nre.s(R, f) - I'I'i-esCR, 01^ are usu- 
ally interpreted as atomic and molecular condensate densities, 
respectively. 



2. Two component model Hamiltonian 

Equations j99> and JlOOI l can be derived formally as equa- 
tions of motion for the mean fields 



4'(x,0 = <'A(x))r 
>I'res(R,f) = <'Ares(R)>r 



(103) 
(104) 



in the Hartree approximation. The corresponding two compo- 
nent model Hamiltonian is given by: 



H 



1 



dx (A"'"(x) 



2m 

+ J ^/r<aL(R) 

+ J dR^rUR) 



2 |„|2 



-4;;^^^---i«i 



dEy, 



dB 



■(fires) 



«A(x) 

<Ares(R) 
[fi(f) - fio]'Ares(R) 



dx (/r"''(x)(/r' (x)l/r(x)i/r(x) 



+ ^fes J fl!x [(AL(x)<A(x)<A(x) + h.c] . (105) 

Here the field operators i/'(x) and i/'ies(R) fulfil the usual 
Bose commutation relations and commute among the differ- 
ent species: 



[lA(x),lAres(R)] =0 
[^(X),^L(R)] =0. 

3. Deficits of the two level mean field approach 



(106) 
(107) 



Despite its common use in studies of molecular association 
in dilute Bose-Einstein condensates, the two level mean field 
model in Eqs. \99\ and ( I100> and its common interpretation 
are subject to two serious deficits. First, in the derivation of 
the mean field equations we have shown that the mean field 
*i',es(R, t) is associated with a pair of atoms in the closed chan- 
nel resonance state 0res('") [see Eq. (I97ll1. Figure|3lclearly re- 
veals that, in general, such an atom pair cannot be associated 
with a molecular bound state because in the close vicinity of 
the 100 mT Feshbach resonance of '^^Rb, as well as asymp- 
totically far from it, the wave function of the highest excited 
vibrational bound state is dominated by its component in the 
open channel |62]. 

The interpretation of Mes(0 = / '^R l*I'res(R, in terms of 
the population in the resonance state (as opposed to the num- 
ber of bound molecules) is a direct consequence of the com- 
mutation relation in Eq. Mifll and of the form of the Hamil- 
tonian in Eq. ( I105> . rather than an artifact of the derivation of 



the mean field equations \99\ and ( I100> . In fact, according to 
Eq. jl07> . a localised pair of atoms created from the vacuum 
state |vac) by the field operator ip\g^(R) is always orthogonal 
to the localised state of any atom pair created by i^' (y)i/r' (x) 
in the open channel; i.e. 



<vac|iA(x)(A(y)iAL(R)|vac) = 0. 



(108) 



The two states of the atom pairs, therefore, correspond to dif- 
ferent asymptotic scattering channels. It can also be verified 
from the Hamiltonian in Eq. MQ51 that any physical two-body 
state associated with i/'i'es(R)|vac) is not stationary with respect 
to the time evolution in the relative motion of the atom pairs. 
The off diagonal coupling in Eq. ( I105l l leads to a decay into 
the open channel, which is determined by the coupling con- 
stant gres- The energy £'res(fi) = [%^(fires)] (fi - fio), which 

can be associated with (/'res(R)lvac), is linear in the magnetic 
field strength and corresponds to the resonance state rather 
than to the bound molecular states (cf. Fig.O. The off diago- 
nal coupling constant gres then determines the corresponding 
decay width. 

Although it seems conceivable to measure the closed chan- 
nel resonance state population M-es(0 - J dR |*l'res(R, 01^^ 
several present day experiments on molecular association in 
Bose-Einstein condensates have clearly revealed the multi- 
channel nature of the bound states. The experiments in 
Refs. |'7,"8'l, for instance, have determined the near resonant 
universal form £'b(fi) = -fi^/ ^m[a(B)]^^ of the binding energy 
in Eq. j39L corresponding to a bound state wave function of 
the universal form of Eq. ( I40> . which is dominated by its com- 
ponent in the asymptotic open channel (cf., also, Fig.|^. The 
experiment in Ref. 1 10] has determined the change of the mag- 
netic moment of the molecular bound states in dependence on 
the magnetic field strength B, which is a consequence of the 
nonlinear dependence of the binding energy on B (see Fig.|5Jl. 
The magnetic moment of the resonance state, however, is in- 
dependent of B. These experiments clearly reveal that the 
number of bound molecules in Eq. (|^} is the relevant physical 
quantity rather than the resonance state population. 

The second major deficit of the mean field equations j99t 
and ( I100> in the description of the adiabatic association tech- 
nique is related to their two level nature; the only configura- 
tions of two atoms are a pair of condensate atoms and the res- 
onance state. Neither the near resonant diatomic bound states 
nor the continuum states can be represented solely in terms of 
these two configurations. In particular, any coupling between 
the initial atomic condensate and the quasi continuum of states 
above the dissociation threshold energy of the open channel 
is ruled out. The two level nature of Eqs. j99t and jlOO> is 
a consequence of the application of the Markov approxima- 
tion to Eqs. iE9i and ( I90l l. This indicates that the presuppo- 
sitions leading to the Markov approximation are violated dur- 
ing the passage across a Feshbach resonance. Consequently, 
Eqs. j99t and JlOOI l fail quantitatively and even qualitatively 
in the dynamic description of the near resonant atom-molecule 
coherence experiments in Refs. |7, 8]. The most striking con- 
sequence of their two level nature, however, consists in the 
insensitivity of Eqs. ( I99> and dlOOt with respect to the ramp 
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direction. The approach thus predicts molecular production 
even for an upward ramp across the 100 mT Feshbach reso- 
nance of '^^Rb, although there is no energetically accessible 
diatomic vibrational bound state on the high field side of the 
resonance (cf. Fig. |2j- The reason for this failure is the ab- 
sence of continuum states that the closed channel resonance 
state could decay into. The results of Subsection lll El suggest. 
however, that the two level mean field approach may recover 
the far resonant asymptotic molecular population in a down- 
ward ramp of the magnetic field strength across the 100 mT 
Feshbach resonance of **^Rb. 



4. Applications of the two component model Hamiltonian beyond 
mean field 

When applied beyond mean 

field 00 El, 

the model 

Hamiltonian in Eq. ( I105> leads to ultraviolet singularities that 
are related to the energy independence of the coupling con- 
stants gbg and gres, i-C- the associated contact potentials cru- 
cially lack any spatial extent. The approach in Ref. |631 cir- 
cumvents this ultraviolet problem in terms of an energy cutoff, 
which is adjusted in such a way that the Hamiltonian, when 
applied to only two atoms, recovers the exact binding energy 
of the highest excited vibrational bound state. This beyond 
mean field approach has been applied recently to the exper- 
iments in Ref. \T^. We expect this approach to give results 
similar to those of the first order micros copic quan tum^ d y - 
namics approach. It has been pointed out [45, 53, 66, 
however, that the interpretation in Ref. 1 63 j of the closed chan- 
nel resonance state population A^resCO in terms of the number 
of bound molecules (i.e. the number of undetected atoms in 
Ref. 1 7]) is not applicable. 

Some approaches |67||, l68ll have suggested curing this 
deficit, by multiplying the closed channel resonance state pop- 
ulation Nresit), as determined in Ref. |63], with a magnetic 
field dependent factor, termed the wave function renormalisa- 
tion factor, which accounts for the component of the bound 
state wave function in the asymptotic open scattering chan- 
nel. A short calculation using Eq. ( I30> shows that the wave 
function renormalisation factor is equivalent to the squared 
normalisation factor 
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abg m[a(B)Y 



(109) 



of the two channel bound state wave function in Eq. ( I29> . pro- 
vided that the magnetic field strength is sufficiently close to 
the resonance position that the bound state wave function and 
its binding energy are universal (see Subsection lll Ct . A com- 
parison between the closed channel resonance state popula- 
tion in Ref. 1 63 ] , revised by multiplication with the right hand 
side of Eq. jl09> (as provided by Ref. l'67'l), indeed, largely 
recovers the magnitude of the molecular fractions predicted 
by the first order microscopic quantum dynamics approach in 
Refs. j4^l53ll . The revised resonance state population, how- 
ever, shows unphysical oscillations between the populations 
of bound and free atoms (cf. Fig. 3 in Ref. |63]), for magnetic 
field strengths at which the spatial extent of the highest excited 



vibrational molecular bound state is by far smaller than the 
mean interatomic distance in the dilute gas. This unphysical 
behaviour is due to the fact that not only the bound state wave 
function in Eq. j29> . but also the continuum wave functions 
in Eqs. ( I16t and illi . have a closed channel resonance state 
component. As a consequence, there is no general relation- 
ship between the closed channel resonance state population of 
a gas and the number of bound molecules. 

The general relationship between the number of bound 
molecules and the two-body correlation function in Eq. ( I87t . 
however, is applicable also to the approach in Ref. 1 63 ] . Equa- 
tion (|^) applies the appropriate quantum mechanical observ- 
able, and will therefore automatically determine the mea- 
surable fraction of bound molecules in a gas, provided that 
molecules identified as separate entities are at all a reasonable 
concept [53]. 



C. Comparison between different approaches 

We have shown in Subsection llll AI that, during a linear 
ramp of the magnetic field strength across a Feshbach reso- 
nance, the full microscopic quantum dynamics approach can 
predict a transfer of atoms from the atomic condensate not 
only to the bound, but also to the continuum part of the two- 
body energy spectrum. The transfer into the continuum en- 
ergy states, however, occurred only in upward ramps of the 
magnetic field strength across the 100 mT Feshbach resonance 
of ^^Rb. The two level configuration interaction approach of 
Subsection lll EI restricts the configuration of an atom pair in 
the open channel to the lowest energetic quasi continuum state 
of the background scattering in the atom trap, while in the two 
level mean field model of Subsection llll Bl all atoms in the 
(F = l,mf - +1) hyperfine state are described by a coherent 
mean field. Both two level approaches, therefore, rule out the 
production of correlated pairs of atoms in the quasi continuum 
of two-body energy levels from the outset. 

Equations i62\ and (I63t of the configuration interaction ap- 
proach and the mean field equations ( I99I I and dlOOI l describe 
the coupling between the atomic condensate and the closed 
channel resonance state in essentially the same way, except 
that the phases of the off diagonal coupling terms are different. 
The spatial configuration of the atomic condensate in the con- 
figuration interaction approach, however, is static, while the 
two level mean field approach allows the trapped atomic con- 
densate to access coherent collective excitations with a high 
occupancy of the energy modes. Consequently, the configura- 
tion interaction approach can, at least to some extent, be inter- 
preted as a local density approximation to the two level mean 
field model. In accordance with the derivations in Subsection 
nil Bl both two level approaches can, therefore, be considered 
as the Markov approximation to the first order microscopic 
quantum dynamics approach of Subsection llll AI 

Our previous considerations have shown that there exists 
a qualitative agreement between the approaches, at differ- 
ent levels of approximation, with respect to the prediction of 
molecular formation in linear downward ramps of the mag- 
netic field strength across the 100 mT Feshbach resonance 
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of ^^Rb. The description of the microscopic binary collision 
physics as well as their treatment of the coherent nature of 
the Bose-Einstein condensate, however, differ considerably 
among these approaches. In the following, we shall provide 
quantitative comparisons between their predictions with re- 
spect to the molecular production efficiency of the adiabatic 
association technique. We will study in detail the dependence 
of the magnitude of the molecular fraction on experimentally 
accessible parameters. 



1. Universal properties of the molecular production efficiency of a 
linear ramp of the magnetic field strength 

In the adiabatic association technique with linear ramps of 
the magnetic field strength, the ramp speed ^ controls the in- 
teratomic interactions. Furthermore, the initial number of 
atoms can be varied over a wide range in present experiments, 
and the trap frequencies determine the confining potential of 
a harmonic atom trap. The atom trap, the number of atoms 
as well as their binary interactions control the spectrum of co- 
herent collective excitations of a Bose-Einstein condensate. In 
the following, we shall consider spherically symmetric atom 
traps with a radial trap frequency that we denote by Vho. 

Figure [21 shows the predicted dependence of the remnant 
atomic condensate fraction and of the fraction of atoms as- 
sociated to diatomic molecules in the highest excited vibra- 
tional bound state on the inverse speed of linear downward 
ramps of the magnetic field strength across the 100 mT Fes- 
hbach resonance of **^Rb. We have chosen trap frequencies 
of 100 Hz (a) and of 10 Hz (b) that correspond to a compar- 
atively strongly and a weakly confining harmonic atom trap, 
respectively. The initial number of Ndti) - 50000 condensate 
atoms is fixed in all calculations. The predictions correspond 
to the full first order microscopic quantum dynamics approach 
(MQDA), the nonlinear configuration interaction equations 
(I62> and (I63> (CI), the two level mean field equations ( I99> 
and ( I100> (2 level), and the Landau-Zener asymptotic popula- 
tions in Eqs. ( I59> and ( I60> (LZ). The comparison reveals that 
the approaches predict similar conversion efficiencies with re- 
spect to the association of molecules as well as similar rem- 
nant condensate fractions. Only the saturation behaviour of 
the exponential Landau-Zener (LZ) curves [cf. Eqs. ( 15 9> and 
(I60H diff'ers significantly from all other approaches because it 
corresponds to the asymptotic populations of the hnear equa- 
tions ( I55> and ( I56> . These results suggest a remarkable in- 
sensitivity of the molecular production in linear ramps of the 
magnetic field strength with respect to both the details of the 
microscopic binary collision physics and the coherent nature 
of the Bose-Einstein condensate. 



2. Dependence of the molecular production efficiency on 
experimentally accessible parameters 

We shall now identify the universal physical parameter that 
controls the molecular production efficiency of the adiabatic 
association technique. As shown in Subsection lll El the linear 
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FIG. 12: Predicted remnant condensate fraction and number of atoms 
associated to diatomic molecules in the highest excited vibrational 
bound state after a linear downward ramp of the magnetic field 
strength across the 100 mT Feshbach resonance of '^^Rb in depen- 
dence on the inverse ramp speed. The initial state of the gas corre- 
sponds to a dilute zero temperature Bose-Einstein condensate with 
a scattering length of Obs = 100 OBohr and A' = 50000 atoms in 
spherically symmetric atom traps with the frequencies 100 Hz (a) 
and 10 Hz (b). The abbreviations in the legend indicate the different 
approaches applied in the calculations. These approaches contain 
the first order microscopic quantum dynamics approach (MQDA), 
the nonlinear configuration interaction equations <62t and <63> (CI), 
the two level mean field equations <991 and <100> (2 level), and the 
Landau-Zener asymptotic populations in Eqs. <59t and <60t (LZ). 
The small deviations in the calculated MQDA data of the remnant 
atomic condensate (filled circles) from an entirely smooth curve are 
discussed in Fig. llll 



Rabi flopping model in Eqs. j55> and j56t leads to an asymp- 
totic fraction of diatomic molecules that can be described by a 
single parameter, the Landau-Zener coefficient 5lz- We shall 
show in the following that also in the nonlinear two level ap- 
proaches the asymptotic molecular population is characterised 
by 5lz- To this end, we take Eqs. j99> and jlOOt of the two 
level mean field approach, where it is revealing to introduce a 
natural time scale fnat - (^^)/(^)' well as the harmonic 
oscillator length scale /ho = y/W/(m2nvho). Applying these 
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FIG. 13: Predicted proportion of remnant condensate atoms and pro- 
portion of atoms associated to diatomic molecules in the highest ex- 
cited vibrational bound state after a linear downward ramp of the 
magnetic field strength across the 100 mT Feshbach resonance of 
*^Rb in dependence on the inverse ramp speed. The atom number and 
the angular frequency a)h„ of the atom trap are adjusted to keep the 
nonlinearity parameter k\,g and the quantity aihotmt constant, which 
also fixes the Landau-Zener parameter Si,z- The initial states of the 
gas correspond to a dilute zero temperature Bose-Einstein conden- 
sate with a scattering length of a(,B = 100 OBohr in spherically sym- 
metric atom traps. In each case the shape of the initial condensate 
mode is, in units of 4^, identical to that of a dilute Bose-Einstein 
condensate with 50000 atoms in a spherically symmetric atom trap 
with a frequency of Vho = 10 Hz. The abbreviations in the legend 
indicate the different approaches applied in the calculations. These 
approaches contain the first order microscopic quantum dynamics 
approach (MQDA), and the two level mean field equations <99> and 
fTO (2 level). 



natural time and length scales to the two level mean field equa- 
tions reveals that they can be characterised purely in terms 
of the nonlinearity parameter k^g - A^flbg/4o of the Gross- 
Pitaevskii equation |47], in addition to the quantities Whofnat 
and 



^eff = 



dEr, 



dB 



- (fires) 



{ABy/(dB/dt). 



(110) 



Here is the total number of atoms and Who = ^^Vho is the 
angular frequency of the spherically symmetric atom trap. We 
note that all three of these parameters can in principle be in- 
dependently varied by manipulating the trap frequency, the 
ramp speed, and the total number of atoms. Of these three 
dimensionless quantities, the Landau-Zener parameter 6lz in 
Eq. ( 16 H can be written as a function of just the nonlinearity 
parameter k\,g and Whofnat- Indeed, as we observe in Fig. 1131 
keeping k\,g and Whofnat constant while varying ^eff reveals a 
very weak dependence of the molecular conversion efficiency 
on this variation. Using identical input parameters for the dy- 
namic equation ( I94l l of the first order microscopic quantum 
dynamics approach reveals that this remarkable universality 
of the molecular production during a linear passage across 
a Feshbach resonance is preserved even when the complete 



quasi continuum of excited two-body energy modes is explic- 
itly accounted for by the theory. 



IV. CONCLUSIONS 

We have presented a comprehensive theoretical analysis of 
the adiabatic association to diatomic molecules of initially 
Bose-Einstein condensed **^Rb atoms via magnetic field tun- 
able interactions. In particular, we have considered the sit- 
uation in which a gas of Bose-Einstein condensed atoms in 
the (F - l,mf = -i-l) state is exposed to a homogeneous 
magnetic field whose strength is varied linearly to cross the 
broadest Feshbach resonance at 100 mT, in order to produce 
strongly correlated pairs of atoms in the highest excited vi- 
brational bound state. We have compared the predictions of 
Landau-Zener, nonlinear configuration interaction, and two 
level mean field approaches with the full first order micro- 
scopic quantum dynamics approach, which explicitly includes 
all energy states of two atoms. 

We found that despite decisive differences between these 
many-body approaches with respect to the description of the 
underlying microscopic binary collision dynamics, the pre- 
dicted molecular production efficiencies obtained from linear 
ramps of the magnetic field strength are virtually indepen- 
dent of the level of approximation. We have shown that the 
Landau-Zener coefficient of Eq. ( I61l l is the main parameter 
that controls the molecular production in all theoretical ap- 
proaches. Consequently, the efficiency of molecular produc- 
tion in linear ramps of the magnetic field strength is remark- 
ably insensitive with respect to the details of the binary colli- 
sion dynamics and to the coherent nature of the gas. The adi- 
abatic association of molecules in dilute Bose-Einstein con- 
densed gases is thus subject to universal physical properties 
similar to those that we have identified in the associated two- 
body problem. This indicates that related experiments on the 
formation of molecules as well as their subsequent dissocia- 
tion via linear ramps of the magnetic field strength (see, e.g., 
Ii50i] ') are largely inconclusive with respect to the details of the 
intermediate microscopic binary collision dynamics (cf., also, 
Ref. 154]). 

The universal properties of the molecular production effi- 
ciencies and dissociation spectra reported in this paper are re- 
stricted to linear ramps of the magnetic field strength. The 
interferometric studies of Refs. [2, ^ and their subsequent 
theoretical analysis i45ll63ll64ll . for instance, have clearly re- 
vealed that sequences of linear variations of the magnetic field 
strength lead to molecular production efficiencies beyond the 
predictions of Landau-Zener or two level mean field models. 
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APPENDIX A: LOW ENERGY TWO CHANNEL 
HAMILTONIAN 

In this appendix we shall provide the universal and practi- 
cal model of the low energy two-channel potential matrix that 
we have used to determine the static as well as the dynamic 
properties of the resonance enhanced scattering in SectionsHll 
and|ni] Based on the general form of the energy states in Sub- 
section irr51 we shall determine the explicit parameters of the 
low energy potentials that correspond to the 100 mT Feshbach 
resonance of '^^Rb. 



1. Low energy background scattering potential 

The two-channel bound and continuum energy states in 
Subsection IIIBI depend on the complete Green's function 
Gbg(z) of the background scattering in Eq. (|5Jl. Here z is the 
complex parameter in Eqs. Mil and Eq. \29l . respectively, 
which determines the asymptotic form of the continuum and 
bound state wave functions at large interatomic distances. The 
Green's function Gbg(z) is determined by all bound and con- 
tinuum energy states associated with the background scatter- 
ing. The detailed form of the binary interaction potential is not 
resolved in ultracold collisions (cf. discussion in 111 B 11 1. At 
magnetic field strengths asymptotically far from the position 
of a Feshbach resonance, the low energy scattering observ- 
ables can be determined in terms of just two parameters, the 
background scattering length Obg and the binding energy E^i 
(see Fig. |5J of the highest excited vibrational bound state of 
the background scattering potential Vbg |37, 69]. This univer- 
sality allows us to choose a potential model of Vbg that recov- 
ers the exact scattering length Obg as well as the exact binding 
energy E-\. In the following, we thus use a most convenient 
separable potential energy operator, of the form fTOll 



Vbg = Urbg>fbgCrbgl, 



(Al) 



to determine Gbg(z). We choose the arbitrary but convenient 
Gaussian form of the function 



(PlYbg) = Xhiip) 



exp 



(A2) 



in momentum space. We shall then adjust the amplitude ^bg 
and the range parameter cTbg in such a way that they exactly 
recover Obg as well as (see Fig.|5}. 

The choice of a separable potential energy operator in 
Eq. JA1> allows us to determine the bound state 0_i and the 
continuum energy states of the background scattering in an 



analytic form. The bound state 0_i and its binding energy £_! 
are determined by the integral form of the Schrodinger equa- 
tion, which for the separable potential in Eq. (lAH is given by: 

|0-i> = Go(£-i)Vbgl0-i> = Go(£-i)brbg)ftg(;irbgl<^-i>. (A3) 

Here Go(z) = [z - [-fi^'^^lmy^ is the free Green's function, 
i.e. the Green's function of the relative motion of an atom pair 
in the absence of interactions. The factor (;\fbgl0-i) in Eq. ( IA3t 
is determined by the unit normalisation of 0_i . The as yet un- 
determined binding energy E^x can be obtained by multiply- 
ing Eq. JA3> by {xhi\ from the left. This yields: 

l-^bgCtbg|Go(£-i)L;tbg> = 0. (A4) 

A short calculation using the Gaussian form of ;^fbg in Eq. JA2> 
determines the matrix element of the free Green's function in 
Eq. (EU by: 



Crbg|Go(£'_i)[^bg) = 



(A5) 

Here erf(x) = ^ C e " du is the error function with the ar- 

^/tt JO 

gumentx= y/mlE-ilcrhg/ft. 

In addition to Eq. JA4l i. the second condition that is needed 
to relate the two parameters ^bg and ixbg of the separable po- 
tential operator to the background scattering length and the 
binding energy £■_] is provided by the continuum energy wave 
functions, or equivalently, by the T matrix. The T matrix 
of the background scattering is determined by the Lippmann- 
Schwinger equation i3^ : 

Tbgiz) = Vbg + VbgGo(z)rbg(z). (A6) 

For a separable potential energy operator Eq. jA6> can be 
solved by iteration, which yields the sum of the Born series: 

lrbg)^bgCrbgl 



7'bg(z) = VbgXN2)^bJ = r 



■ ^bg(Xbg\Go(z)\xbg}' 



(A7) 



The background scattering length is then determined in terms 
of the T matrix in Eq. iAl\ by 



Clbg 



4;rfi^ fbg 



i-^bgCrbglGo(O)Ufbg) 



(A8) 

Here |0) is the zero momentum plane wave of the relative mo- 
tion of an atom pair The denominator on the right hand side 
of Eq. (|A8) 

can be obtained directly from Eq. iA5\ by replac- 
ing the energy argument E^i by 0. This yields: 



dbg 



m ^ 



(A9) 



l + 4i^^bg/(V^crbg)' 

Equation jA9l l can be used to eliminate ^bg from Eq. <A4l l. 
which, in turn, determines the range parameter crbg in terms 
of the background scattering length Obg and of the binding en- 
ergy £■_!. Given the range parameter crbg the remaining am- 
plitude ^bg can then be obtained from Eq. jA9l l. For the '^^Rb 
parameters Ubg =100 flBohr and E^i = -h x 23 MHz, as ob- 
tained from Fig. ^ this yields crbg = 42.90753599 flBohr and 
OT^bg/ (47rri2) = -317.5649079 flBo 



Bohr- 
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2. Off diagonal coupling 

In order to recover the universal properties of the near res- 
onant bound states in Subsection 111 CI (cf . Fig.O, the transi- 
tion probabilities in Subsection lll DK cf. Fig.|5jl, as well as the 
time evolution operator in Eq. i94l . it is sufficient to use a sin- 
gle channel Hamiltonian, as explained in Subsection lll CI To 
this end, we have extended the potential model in Eq. (lAH 
to recover not only the background scattering, but also the 
scattering length a(B). We have thus adjusted the separable 
potential energy operator in Eq. iAl\ . at each magnetic field 
strength B, to the magnetic field dependent scattering length 
a(B) = flbe [1 - (^B)/(B - Bo)] and to the energy E^i (cf., 
also, Ref. |45|). 

The dissociation spectra in Subsection III Fl have been ob- 
tained from a two channel description. A two channel ap- 
proach recovers not only the exact magnetic field dependence 
of the scattering length, but it also accurately describes the 
binding energies of the multichannel bound states over a much 
wider range of magnetic field strengths (see Fig.|2} than a sin- 
gle channel treatment. Apart from the background scattering, 
a two channel description of resonance enhanced collisions 
requires us to specify the coupling between the open and the 
closed channels. Equations i25i and ( I26> relate the off diago- 
nal coupling element W(r) of the general two channel Hamil- 
tonian in Eq. to the resonance width (AB) as well as to the 
shift Bo - Bres- Although the resonance shift is not directly 
measurable, Eq. iiH relates it to the van der Waals dispersion 
coefficient and to the width (AB), which can usually be 
determined from experimental data. In the pole approxima- 
tion to the closed channel Green's function [cf. Eq. (I15H. the 
closed channel part of the two channel Hamiltonian in Eq. 
reduces to Eq. (I20> . which restricts the state of an atom pair 
in the closed channels to the resonance state 0ies. Given the 
linear dependence of the energy £'res(B), associated with (^res, 
on the magnetic field strength B, the closed channel part of 
the two channel Hamiltonian is determined completely by the 
slope [^^(fires)] of the Feshbach resonance (cf. Subsection 

The only remaining quantity that needs to be determined is 
the product W(r)<presif), which provides the interchannel cou- 
pling in the pole approximation. Similar to the choice of the 
separable potential in Eq. iAl\ . we shall use a two parame- 
ter description of the interchannel coupling, in terms of a real 
amplitude ^ and a range parameter cr, to recover the width and 
the shift of the Feshbach resonance. This leads to the ansatz 



(AlO) 



where we have chosen the arbitrary but convenient Gaussian 
form of the function 



exp 



<plr) = xip) = 



(-*) 



(All) 



in momentum space. In the following, we shall adjust the pa- 



rameters f and cr in such a way that the two channel Hamilto- 
nian recovers the resonance width {AB) and the shift Z?o - Bres 
via Eqs. \25\ and \26\ . respectively. To this end, we shall first 
determine the zero energy continuum state (f}^ ^ in Eq. \25\ and 
the zero energy Green's function Gbg(O) in Eq. \26\ in terms of 
the T matrix associated with the background scattering, which 
is given by Eq. \A6\ . The continuum energy states associated 
with the relative momentum p fulfil the Lippmann-Sch winger 
equation JH (£ = p^lm): 



I0i,^'> = IP> + Go(£ + /0)rbg(£ + /0)|p). 



(A12) 



Here the energy argument "z = £ + /O" ensures that the wave 
function 0p^^(r) has the long range asymptotic form of Eq. 
Furthermore, the Green's function Gbg(z) is completely deter- 
mined by the T matrix in Eq. \A6\ in terms of the general 
formula L38.1 : 



Gbg(z) = G„(z) + Go(z)rbg(z)Go(z). 



(A13) 



In the separable potential approach the exact T matrix of the 
background scattering is given by Eq. jA7l l. so that the subse- 
quent determination of the right hand sides of Eqs. \25\ and 
\26\ involves just the evaluation of matrix elements of the zero 
energy free Green's function Go(0) between the wave func- 
tions ;t'bg and;t'. Given the Gaussian form of these wave func- 
tions in Eqs. JA2l i and \A\ 11 1, the matrix elements can be read- 
ily determined analytically in complete analogy to Eq. \A5\ . 
This yields 



{AB) 



flbg 



for the resonance width and 



(A14) 



Bo - Bres = (AB) 



«bg 



(A15) 



for the shift. Here we have introduced the mean range param- 
eter 



^=V^V + <)- (A16) 

Inserting the experimentally known width of the resonance 
and Eq. ( I32t for the shift into the left hand side of Eqs. ( IA14t 
and jA15> . respectively, then determines ^ and cr in terms of 
(AB), and the slope of the resonance. Using (AB) = 0.02 
mT 144], Ce = 4660 a.u. |42] and [^(Bi-es)] = x 38 
MHz/mT for the 100 mT Feshbach resonance of **^Rb, the 
parameters of the interchannel coupling can be summarised 
as cr = 21.50035463 OBohi- and mi^^ l{\Tih!'(T) = /z x 8.0536126 
MHz. 
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